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STRUCTURAL PETROLOGY OF CAMBRO-ORDOVICIAN 
LIMESTONES OF THE NORTH-WEST HIGHLANDS 
OF SCOTLAND 


C. B. CRAMPTON 


ABSTRACT. Petrofabric analyses of dolomite and calcite marbles in the Durness Lime- 
stone from Loch Broom to Loch Glencoul in the North-West Highlands of Scotland are 
described. The preferred orientation of compression axes (and initial [000] ]-axes) 
constructed for twinned {()22]}-lamellae of dolomite measured in these limestones shows 
a transition from partial girdles characterised by two mutually perpendicular and promi- 
nent maxima to complete girdles disposed vertically either south-east and north-west or 
— and north-east. A tentative kinematic interpretation of these fabrics is at- 
tempted. 
INTRODUCTION 

Both calcite and dolomite have been investigated under conditions of ex- 
perimental deformation (Griggs, Turner, Borg and Sosoka, 1953; Handin and 
Fairbairn, 1953; Turner, Griggs, Heard and Weiss, 1954). Some aspects of 
certain fabrics of naturally deformed calcite and dolomite marbles have been 
interpreted on the basis of the experimental data (Turner, 1953; Gilmour and 
Carman, 1954; Crampton, 1956). Petrofabric analyses of chiefly dolomite 
and some calcite marbles in the Durness Limestone succession of the North- 
West Highlands of Scotland from Loch Broom to Loch Glencoul are here 
described. 

Situation and general features of analyzed specimens.—Loch Broom lime- 
stone: collected from the north-east cliffs of Loch Broom at Gadcaisceig, one 
mile south-east of Ullapool, not far below the Moine Thrust. The bedding 
plane dips 13° south-east. 

Lochan Fasaidh limestone: collected just below the Moine Thrust a quar- 
ter of a mile north of Lochan Fasaidh. The foliation dips 47° south-south-west. 
Dip and bearings of sections cut for fabric analysis are 44° /337°, 61°, 085° 
and 47°/205°. 

Loch Urigill limestone: collected from the shore of Loch Urigill at the 
base of the north-east slopes of Druim na Doire Duibhe, mid-way between the 
thrust planes and sole. The foliation is strongly ruckled, but at the point of 
collection dips 44° west. Dip and bearing of sections cut for fabric analysis 
are 85°/176°, 52°/088° and 35° /280°. 

Elphin limestone: collected from the south-west slopes of Cnoc Breac just 
east of Elphin and the sole of the thrusts. The rudimentary foliation dips 30° 
west-north-west. Dip and bearing of sections cut for fabric analysis are 30° 
293°, 70°/132° and 73°/032°. 

Ledbeg limestone: collected from the banks of the Ledbeg River near 
Ledbeg, close to the sole of the thrusts. The rock is cut by a set of parallel 
calcite veins along which splitting tends to occur thus giving rise to a well 
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defined joint system. The foliation dips 14° east, and carries a fine lineation 
dipping 11° north-east. This lineation lies normal to the calcite veins. Dip and 
bearing of sections cut for fabric analysis are 82°/325°, 81°/220° and 14°, 
032 

Cnoc Chaornaidh limestone: collected from the lower slopes of Cnoc 
Chaornaidh, alongside the road to Benmore Lodge and near the Benmore 
Thrust. The strongly ruckled foliation at the point of collection dips 20° west. 
Dip and bearing of sections cut for fabric analysis are 90° 311°, 82°/135 
and 08° /315 

Loch Glencoul limestone: collected immediately below the Glencoul 
Thrust plane near the shores of Loch Glencoul. The bedding dips 20° south- 
east. 

In thin section the Loch Broom limestone shows angular dolomite and 
quartz grains of sedimentary character. A little nearer the Moine Thrust plane 
there is some evidence of recrystallization, and the dolomite grains are larger 
and flattened normal to the optic axis. Against the lower surface of the Glen- 
coul Thrust plane the Loch Glencoul limestone shows a mosaic of dolomite 
crystals intersected by zones of fracture, which sometimes contain lenticles of 
quartz. 

Some specimens were collected from recrystallized masses within the 
metamorphic aureole of the Cnoc na Sroine syenite. Within the Lochan 
Fasaidh and Elphin limestones. clusters of dolomite crystals are enclosed by 
granulated dolomite. In the Ledbeg limestone individual dolomite porphyro- 
blasts with irregular outline are isolated by the dolomite granules (plate 1), 
while within the Cnoc Chaornaidh limestone large flattened dolomite por- 
phyroblasts are enveloped by granoblastic calcite and brucite (plate 2). 

LOCHAN FASAIDH LIMESTONE 

In thin section clusters of dolomite porphyroblasts are enclosed by granu- 
lar dolomite. A little marginal granulation suggests slight external rotation of 
some porphyroblasts since crystallization. Most grains possess {0221}- 
lamellae and {1011}-cleavages, these planes being identified by their angular 
and zonal relationships. Slight internal distortion of some porphyroblasts is 
indicated by undulose extinction and flexed lamellae and cleavages. 

Petrofabric analysis——Because in no case were nearly all twinned 
lamellae in the fabric accessible for measurement in one section only, three 
mutually perpendicular sections were cut from all specimens described within 
this study and used for investigation of preferred orientation. 200 grains were 
measured in each section. Unless the figures have been presented in the 
horizontal plane showing peripheral geographical co-ordinates, the plane of 
projection of each diagram coincides with the plane of that section yielding 
the highest proportion of twinned grains. Plotted points representing measure- 
ments made in the other two sections have been rotated to the plane of the 
former. 

Within this particular specimen only two sections revealed a significant 
proportion of twinned grains. Figure 1 shows the distribution of 400 twinned 
and untwinned grains measured in these two sections. The restricted range of 
orientation of the initial [0001 ]-axis of those grains with optically identifiable 
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(0221)}-twins contrasts sharply with apparently random distribution of [0001] 
of twinned and untwinned grains considered collectively. Hence the range of 
lattice orientation does not limit the range of possible orientation of {0221}- 
twins within the fabric. 


Fig. 1. Loch Fasaidh limestone. Dolomite: 75 initial [000]]-axes of twinned grains 
(large circles), and 325 [00]]-axes of untwinned grains (small circles), measured in two 
mutually perpendicular sections. Foliation left to right; E, S indicate compass points. All 
figures constructed in lower hemisphere of equal-area projection. 

Fig. 2. Loch Fasaidh limestone. Dolomite: compression axes (63 from the section 
parallel to the plane of projectidn of the diagram (small circles), and 18 rotated from 
the other section (large circles), and tension axes (crosses) constructed for 81 twinned 
{0221 )}-lamellae measured in two mutually perpendicular sections. Foliation left to right; 


E, S indicate conypass points. 


Most grains showing twinning possess only one set of twin lamellae. The 
preferred orientation of {0221}-poles for twinned lamellae define two sta- 
tistical s-planes diverging 45° either side of the foliation. During experimental 
deformation of dolomite (Turner, Griggs, Heard and Weiss, 1954) twin 
gliding along {0221}-lamellae was restricted to those grains in which [0001] 
was originally inclined to the axis of compression at a small angle, the normal 
to the intersection of {0221} with {0001} being the glide direction and the 
sense of shear such that upper layers of the lattice are displaced away from 
the initial [000] ]-axis. A compression and a tension axis most favorably 
oriented to effect twin gliding on {0221} were constructed graphically for all 
twinned lamellae measured in both of the mutually perpendicular sections. 
Figure 2 shows the resultant distribution. 

The preferred orientation of compression axes follows closely that of 
initial [0001 ]-axes of twinned grains because of twinning geometry. Compres- 
sion axes constructed from the lower proportion of measured grains with 
(0221}-twins in the section normal to the plane of projection of the diagram 
confirm the maximal concentration of such axes normal to the foliation. The 
collective preferred orientation of tension axes is, however, slightly weaker. 

Interpretation of fabric data.—Because of the slightly higher preferred 
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orientation of compression axes, more of the observed twinning could have 
been produced by a simple compressive stress operating normal to the folia- 
tion than by a simple tensile stress. The preferred orientation of initial 
[0001 ]-axes of twinned grains within a maximal concentration normal to the 
foliation (fig. 1), on the basis of the experimental data, is consistent with 
these tentative conclusions. Such compression would be directed north-north- 
east and south-south-west, nearly parallel to the line of the thrust outcrop. 
This fabric closely resembles that of the Loch Shin Limestone (Crampton, 
1956) except that the textures and the partial girdle of initial [0001 ]-axes 
suggests minor intergranular rotation in addition to twinning in response to 


stress. 


LOCH URIGILL LIMESTONE 

In thin section the rock consists of granoblastic dolomite dissected by ir- 
regular strings of dolomite porphyroblasts. There is little evidence of inter- 
granular rotation. Undulatory extinction is infrequently present and {0221}- 
lamellae and {1011}-cleavages tend to be straight and widely spaced. Sutural 
contacts between grains are usual. 

Petrofabric analysis.—As in the previous and all following examples 
discussed (with one exception) there is random distribution of [0001] of 
twinned and untwinned grains considered collectively. Hence there is no re- 
striction to possible orientation of twinned lamellae other than orientation of 
the stress system to which twinning may be considered a response. 


Figure 3 shows the distribution within two coincident maxima nearly 


Fig. 4 

Fig. 3. Loch Urigill limestone. Dolomite: compression axes (small circles), and 
tension axes (crosses), constructed for 68 twinned {022]}-lamellae measured in the 
section parallel to the plane of projection of the diagram. Foliation left to right; E, S 
indicate compass points. 

Fig. 4. Loch Urigill limestone. Dolomite: 20 compression axes (large circles) ro- 
tated from one section, and 17 compression axes (small circles) rotated from the other 
section normal to the plane of projection of the diagram, and 37 tension axes (crosses). 
constructed for twinned {022]}-lamellae measured in these two sections. Foliation left 
to right; S, E indicate compass points. 
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parallel and normal to the foliation of compression and tension axes con- 
structed for {0221}-lamellae of twinned grains measured in the section reveal- 
ing the highest proportion of grains with {0221}-twins. Figure 4 shows the 
distribution of these same axes gained from the other two of the three mu- 
tually perpendicular sections cut from each specimen. A comparison of the 
two figures confirms the reality of the two maxima of compression axes at 
right angles to one another, which are coincident with the preferred orienta- 
tion of initial [0001 ]-axes of these same grains. A similar comparison regard- 
ing tension axes reveals that these are distributed within two girdles at right 
angles to one another, each being normal to one of the maxima of compres- 
sion axes. It follows that twinned lamellae poles will tend to be distributed 
along small circles arranged around each of the two maxima of initial [0001 ]- 
axes. 

Interpretation of fabric data.—Because of the stronger and homogeneous 
preferred orientation of compression axes most observed twinning could have 
been produced by compressive stresses, these being of insufficient intensity to 
cause inter-granular rotation. The two maximal concentrations of compression 
axes are directed, one nearly vertically, and the other nearly horizontally east 
and west. 

Regional data.—Analyses of the Elphin, Loch Glencoul and Loch Broom 
limestones reveal fabrics essentially similar to that described above. Figure 5 
shows the distribution of compression and tension axes constructed for 
{0221}-twins measured in three mutually perpendicular sections cut from the 
Elphin limestone. The two marked concentrations of compression axes are 
quite distinct, one nearly vertical and the other nearly horizontal south-east 


Fig. 5. Ebphin limestone. Dolomite: compression axes (contoured), and tension 
axes (crosses), constructed for 154 twinned {022] \-lamellae measured in three mutually 
perpendicular sections. Plane of diagram is the geographic horizon; NSEW indicate 
compass points; foliation shown as full arc. Contours; %4-1-2-3-5% per 1% area. 
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and north-west, both within a south-east and north-west girdle. The lower 
preferred orientation of tension axes is also apparent, these being distributed 
within two girdles, one nearly horizontal and the other nearly vertical, with 
a maximal concentration at their intersection. 

The maximal concentration of compression axes nearly normal to the 
foliation is common to all these limestone fabrics, but the other concentration 
and also the orientation of the girdle of such axes vary approximately parallel 
or normal to the general trend of the thrust outcrop. Analysis of the Loch 
Broom limestone near the Moine Thrust plane reveals the two maximal con- 
centrations of compression axes within a partial girdle disposed north-north- 
east and south-south-west, contrary to that of the Elphin limestone near the 
sole. The fabric of the Loch Broom limestone shows an exceptional weak 
preferred orientation of [0001] of twinned and untwinned grains considered 
collectively normal to the foliation, and coincident with one of the maxima 
within the much higher preferred orientation of initial [0001 ]-axes of twinned 
grains. Experimental deformation of dolomite (Turner, Griggs, Heard and 
Weiss, 1954) did not produce fabric changes of sufficient magnitude to allow 
definite conclusions, but a slight migration of [0001] towards the axis of 
compression was reported. 


LEDBEG LIMESTONE 

Dolomite is present in two forms. The granoblastic dolomite is segregated 
into alternating layers of varying grain size thus giving rise to the foliation. 
These layers usually maintain a constant orientation except in the region of 
porphyroblasts where they are strongly flexed. The delomite porphyroblasts 
display strongly developed undulatory extinction, and are frequently fractured 
and granulated around their margins. They contain irregularly shaped 
{0221}-lamellae and {1011}-cleavages which are often flexed, the region of 
flexing being characterised by a greater development of twinned lamellae. 
There is sometimes a marked loss of parallelism of {0221}-lamellae within 
each set. These textures are illustrated by plate 1, and suggest that the por- 
phyroblasts have undergone severe inter- and intra-granular rotation since 
their formation. Although slight elongation of the porphyroblasts helps pro- 
duce the lineation, this is mainly defined by elongated brucite lenticles. 

The calcite veins lying normal to the lineation tend to possess sharp 
boundaries with the remainder of the rock. In places these veins are flexed, 
this being accompanied by some granulation of the calcite. Undulose extinc- 
tion is present in some grains. While many {0112}-lamellae are fine, straight 
and widely spaced, some grains possess zones of flexure in which there is a 
profuse development of broad lamellae with diffuse margins. These textures 
suggest limited crystal rotation, this never being as severe as that within the 
dolomite porphyroblasts. 

Analysis of the fabric of dolomite——Figure 6 shows the distribution of 
600 [0001 ]-axes of twinned and untwinned grains measured in three mutually 
perpendicular sections, the plane of projection of the diagram lying normal 
to the lineation and parallel to the calcite veins. If only the initial [0001 ]-axis 
of all measured grains with {0221}-twins is selected there is an obvious pre- 
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PLATE 1 


Ledbeg limestone; isolated and fractured porphyroblast showing flexed and disturbed 
twinned {()22] }-lamellae and {10]])}-cleavages. Polaroids crossed; X 40. 


ferred orientation within a girdle normal to the lineation. The full range of 
lattice orientation does not restrict the range of possible orientation of {0221)}- 
twins within the fabric. 

Despite the irregular outline and sometimes diffuse boundary between 
lamellae, in the majority of porphyroblasts the initial lattice could be dis- 
tinguished from the twinned lattice. Figure 7 shows the distribution of com- 
pression and tension axes constructed from twinned lamellae of grains meas- 
ured in the section normal to the lineation. Both axes apparently lie within 
a girdle normal to the lineation. The girdle of compression axes must follow 
from the distribution of initial [0001 ]-axes of twinned grains shown in figure 
6, but the preferred orientation of tension axes needs further investigation. 

Compression and tension axes determined for twinned lamellae of grains 
measured within the remaining two sections, one parallel to the foliation and 
the other normal to the foliation and veins, have been rotated to the plane of 
projection of figure 7. If the distribution of tension axes within figures 7 and 


8 are considered collectively it can be seen that such axes possess random 
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orientation. Any apparently preferred orientation results from consideration 
of measurements made within one section only. Compression axes shown in 
figures 7 and 8 fall within coincident girdles. It follows geometrically that 
twinned lamellae poles are distributed within a weak girdle normal to the 
lineation. 


Fig. 6 Fig. 7 


Fig. 6. Ledbeg limestone. Dolomite: 150 initial [0001]-axes of twinned grains (large 
circles), and 450 [0001]-axes of untwinned grains (small circles), measured in three mu- 
tually perpendicular sections. Foliation left to right; lineation (L); E, N indicate com- 
pass points. 


Fig. 7. Ledbeg limestone. Dolomite: compression axes (small circles), and ten- 
sion axes (crosses), constructed for 93 twinned {()22] }-lamellae, measured in the section 


parallel to the plane of projection of the diagram. Foliation left to right; lineation (L) ; 
E, N indicate compass points. 


Interpretation of fabric data.—Within figure 12 all measured axes from 
the three mutually perpendicular sections have been rotated to their geo- 
graphical orientation. The contoured distribution of compression axes within 
a vertical girdle disposed south-east and north-west, parallel to the calcite 
veins and normal to the foliation and lineation, contrasts with random orienta- 
tion of tension axes. Within the Elphin limestone textures suggesting crystal 
rotation are weakly developed as is the girdle of compression axes shown in 
figure 5. In the Ledbeg limestone the textures of dolomite suggest that intense 
crystal rotation has occurred subsequent to crystallisation and much of the 
twinning. Although the two maximal concentrations of compression axes in 
figure 12 retain a similar orientation compared with figure 5, the well de- 
veloped girdle is now the prominent feature of the preferred orientation. 

Interpretation of fabric of calcite——Only two sections revealed a signifi- 
cant proportion of twinned grains within the total 200 grains measured in 
each section. Figure 9 shows the distribution of compression and tension axes 
constructed from {0112}-lamellae of these twinned grains. Within the girdle 
of compression axes normal to the lineation there is a weighted maximum near 
the foliation. The girdle of tension axes is less well developed. 


St 
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Fig. 8. Ledbeg limestone. Dolomite: 40 compression axes (small circles) rotated 
from one section, and 29 compression axes (large circles) rotated from the other section 
normal to the plane of projection of the diagram, and 69 tension axes (crosses), con- 
structed for twinned {()22]}-lamellae measured in these two sections. Foliation left to 
right; lineation (L); E, N indicate compass points. 

Fig. 9. Ledbeg limestone. Calcite: compression axes (small circles), and tension 
axes (crosses), constructed for 120 twinned {Q]12}-lamellae measured in two of the 
three n-utually perpendicular sections. 30 acute glide line bisectors (large circles) con- 
structed for twinned lamellae of doublets superposed. Foliation left to right; lineation 
(L); E, N indicate compass points. 


Of 90 twinned grains, 30 are doublets. This proportion of doublets need 
not imply that tension was a significant stress during deformation since 
limited crystal rotation suggested by the textures would probably bring more 
than one {0112}-iamella in a grain into a position conducive to twinning. The 
orientation of acute bisectors of glide lines constructed from the {0112}- 
lamellae of these doublets is superposed upon figure 9. When the development 
of twinning has been caused by a compressive stress, this should be approxi- 
mately parallel to the acute bisector of the two sets of lamellae in those grains 
that show double twinning (Turner, 1953). Most of the acute glide line bi- 
sectors shown in figure 9 are located within the broad maximum of compres- 
sion axes. 

Comparison of dolomite and calcite fabrics—Symmetrologically the 
fabrics of dolomite and calcite are similar. The parallel veins containing 
calcite cut the dolomite porphyroblasts, and because they bear a suitable rela- 
tionship with the lineation and preferred orientation of compression axes, are 
interpreted as ac joints which formed after the plastic deformation of dolomite 
had ceased. Two deformation phases must therefore be postulated. Since 
dolomite is considerably stronger than calcite and the textures suggesting 
crystal rotation show weaker development within calcite, the earlier deforma- 
tion phase characterised by plastic deformation of dolomite was far more 
severe than the late minor deformation phase during which calcite experienced 
plastic deformation, after the formation of ac joints. Twinning within calcite 
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could have been produced by a compressive stress directed south-east and 
north-west. 


CNOC CHAORNAIDH LIMESTONE 


An extremely rough appearance on weathered surfaces arises from 
dolomite porphyroblasts standing out from enclosing calcite. In thin section 
brucite forms fibrous aggregates displaying strong undulatory extinction. The 
isolated porphyroblasts are flattened normal to the optic axis, and show those 
textures previously described suggesting limited crystal rotation. The boun- 
daries between adjacent grains within the calcite mosaic are not always dis- 
tinct. These grains are chalky in appearance and {0112}-lamellae are pro- 
fusely developed, such textures having been described by Knopf (1949) and 
Turner and Ch’ih (1951) in calcite experimentally deformed at low tempera- 
tures. The majority of calcite grains are doublets, some of which show one set 


PLATE 2 


Cnoc Chaornaidh limestone; isolated and flattened porphyroblasts showing twinned 
\-lamellae, enveloped by calcite and brucite, Polaroids crossed; 40. 
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profusely developed while the other set consists of widely spaced, relatively 
very thick lamellae with crenulate margins. 

Interpretation of fabric of dolomite-—As in previous examples, [0001] 
of 600 grains measured in three mutually perpendicular sections shows ran- 
dom orientation, while selection of twinned grains reveals a coincident pre- 
ferred orientation of the initial [0001 ]-axis and compression axis. Compres- 
sion and tension axes constructed from twinned {0221}-lamellae measured 
in all three sections have been rotated to their geographical orientation and 
the distribution of compression axes contoured in figure 11. As in the case of 
the dolomite fabric of the Ledbeg limestone, in contrast to random orientation 
of tension axes, compression axes are distributed within a vertical girdle, in 
this case with the contrary orientation south-west and north-east. Closer in- 
spection of this preferred orientation reveals crossed girdles with slight 
divergence, a feature only partially developed within the Ledbeg limestone 
(figure 12). 

At the intersection of the two girdles, the vertical maximal concentration 
of compression axes normal to the foliation, characteristic of most dolomite 
fabrics, is still retained. The horizontally disposed maximal concentration 
within the girdle is, however, replaced by two sub-maxima, one in each of the 
crossed girdles. 


Fig. 10. Cnoc Chaornaidh limestone. Calcite: 66 acute glide line bisectors (small 
circles), and 132 tension axes (crosses), constructed for twinned {0112} -lamellae of 66 
doublets within a total 100 measured twinned grains located within the section parallel to 
the plane of projection of the diagram. Foliation left to right; N, W indicate compass 
points. 

Fig. 11. Cnoec Chaornaidh limestone. Dolomite: compression axes (contoured), and 
tension axes (crosses), constructed for 223 twinned {()22] }-lamellae measured with- 
in three mutually perpendicular sections. Plane of diagram is the geographic horizon; 
NSEW indicate compass points; foliation shown as full are. Contours; 44-1-2-45% per 
1% area. 
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Textures suggesting crystal rotation are more strongly developed within 
dolomite of the Ledbeg limestone. Presumably rotation of the dolomite por- 
phyroblasts would be more easily accomplished in the softer calcite of the 
Cnoc Chaordaidh limestone than within the granular dolomite of the Ledbeg 
limestone. 

Interpretation of the fabric of calcite—Two sections reveal-d very few 
twinned grains. Within the section parallel to the plane of projection of the 
diagram all grains have twinned. Thus calcite had recrystallized with a pre- 
ferred orientation of [0001]-axes severely restricting the range of possible 
orientation of {0112}-twins within the fabric. The distribution of acute glide 
fine bisectors and tension axes constructed from twinned {0112}-lamellae of 
66 doublets of a total 100 measured grains is shown in figure 10. The high 
proportion of doublets implies that a tensile stress could have given rise to 
most of the observed twinning, further confirmation of this being given by the 
lower preferred orientation of acute glide line bisectors. 

Because dolomite porphyroblasts have experienced some crystal rotation 
whereas calcite which is more easily deformed has not, the dolomite fabric 
reflects an earlier and more severe deformation phase. The weighted maxi- 
mum of tension axes of calcite twinned grains coincides with the unoccupied 
area separating the acute angle between the crossed girdles of compression 
axes of dolomite twinned grains. 

Singlet-doublet relationships of dolomite——The geometry of the dolomite 
lattice is such that simple compression applied sub-parallel to [0001] can 
develop high resolved shear stress in the direction and sense suitable for 
twinning on two or all three {0221}-lamellae. Simple compression is con- 
sidered to have been responsible for most observed twinning in the Lochan 
Fasaidh limestone, but only 8 percent of all measured twinned grains were 
doublets, the remainder being singlets. A similar low proportion of doublets 
has been described within the Loch Shin Limestone (Crampton, 1956), and 
occurs in the other limestone fabrics discussed here; for example the Loch 
Urigill limestone 10 percent, Elphin limestone 7 percent, and Ledbeg lime- 
stone 8 percent. The highest proportion of doublets occurs within the fabric 
of the Cnoc Chaornaidh limestone, 24 percent. 

Within this and most other specimens {0221}-lamellae of twinned grains 
are well developed and their twinned character can be determined up to quite 
high angles of tilt thus minimizing the possibility of their being mistakenly 
interpreted as untwinned lamellae. Under comparable conditions of experi- 
mental deformation dolomite was much stronger than calcite and twinning 
far less easily induced (Turner, Griggs, Heard and Weiss, 1954). The absence 
of extensive doublet formation in dolomite might thus be a reflection of this 
greater resistance to plastic deformation by twin gliding, and the wider 
separation of {0221}-lamellae within the dolomite lattice. 


SUMMARY AND CONCLUSIONS 


The calcite fabrics of the Ledbeg and Cnoc Chaornaidh limestones are 
symmetrically similar to the dolomite fabrics within those specimens but be- 
long to a late minor deformation phase. It is of interest that the textures of the 
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Cnoc Chaordaidh limestone suggest cool conditions during deformation by 
comparison with experimental work, and those of the Ledbeg limestone late 
introduction or segregation into veins and subsequent twinning. 


If compression axes of all dolomite fabrics analyzed within this regional 
study are rotated to their geographical orientation and combined in figure 13, 
it is apparent that such axes tend to lie in two vertical girdles, one directed 


S 
Fig. 13 
Fig. 12. Ledbeg limestone. Dolomite: compression axes (contoured), and ten- 
sion axes (crosses), constructed for 162 twinned {()22])}-lamellae measured in three 
mutually perpendicular sections. Plane of diagram is the geographic horizon; NSEW 
indicate compass points; foliation shown as full arc, Contours; 44-1-2-4-55% per 1% area. 
Foliation (S/), calcite veins (S2), and lineation (L) superposed. 


Fig. 13. Composite diagram. Dolomite: 1037 compression axes. Plane of diagram 
is the geographic horizon; NSEW indicate compass points, Contours; 1-2-3% per 1% 
area 


south-east and north-west and the other south-west and north-east, with a 
maximal concentration at the point of intersection of the two girdles. 

A maximal concentration of compression axes (and initial [0001 ]-axes) 
of twinned grains oblique but more usually nearly normal to the foliation is 
present in all the dolomite fabrics. In all cases but one the foliation plane is 
gently inclined to the horizontal and hence this concentration is disposed 
nearly vertically. The exception, the fabric of the Lochan Fasaidh limestone, 
shows a single maximum normal to the foliation (fig. 2), The remaining 
fabrics show a transition from partial girdles characterised by two mutually 
perpendicular and pronounced maxima to complete girdles in which the two 
maxima, although present, are less dominant (figs. 3, 5 and 12). Figures 11 
and 12 show the tendency towards the development of crossed girdles, with 
the maximal concentration at the girdle intersection and the unoccupied area 
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indicating the acute bisectrix of the two girdles near the horizontal plane. 
Figure 11, representing the end of this transition series, shows the reduction 
of one of the two usually typical maxima. 

Attempting a tentative kinematic interpretation of fabric symmetry, if 
the axis of symmetry of the girdles of compression axes or initial [0001 ]-axes 
(the lineation in the fabric of the Ledbeg limestone) is equivalent to the axis 
of external rotation, the deformation (ac) plane suggested by the plane of 
the girdle can be disposed parallel or normal to the general trend of thrust 
outcrops (fig. 13). Thus within the Elphin and Legbed limestones for example 
(figs. 5 and 12) near the sole of the thrusts the plane of rotation would be 
disposed south-east and north-west, which is consistent with the generally ac- 
cepted direction of thrust movements, although only the final stage of de- 
formation is recorded in twin lamellae. In certain places near the thrust 
planes, for example within the Cnoc Chaordaidh (fig. 11) and Loch Broom 
limestones, the plane of rotation would be disposed south-west and north-east, 
suggesting lateral movement at these points. 
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DYNAMIC INTERPRETATION OF THE FABRIC OF 
A DOLOMITE FROM THE MOINE THRUST-ZONE 
IN NORTH-WEST SCOTLAND 
JOHN M. CHRISTIE 


ABSTRACT. Petrofabric studies on naturally and experimentally deformed dolomite 
indicate that the main mechanisms of plastic deformation of dolomite, under most con- 
ditions, are twin-gliding on {0221} and translation-gliding on {0001}. The fabric of a 
dolomite with mylonitic textures is analysed and interpreted in dynamic terms. The grains 
in the rock contain numerous {0221} twin lamellae and also internally rotated lamellae 


of the type observed in experimentally deformed dolomite and designated L,. Compression 
and tension axes inferred from twinned {022]} lamellae are statistically parallel to those 


inferred from L, lamellae. The strain calculated from the amount of rotation of L. 
lamellae varies from 5 to 18 percent for individual grains; it is concluded that these 
figures are probably less than the total post-crystalline strain in the rock. 


INTRODUCTION 
Recent fabric studies of experimentally deformed dolomite rock (Turner, 
Griggs, Heard, and Weiss, 1954; Handin and Fairbairn, 1955) and of single 
crystals of dolomite (Higgs and Handin, 1954) have elucidated the mechan- 
isms of plastic deformation of dolomite at temperatures below 400°C. The 
present paper describes the fabric of a dolomite tectonite from the Moine 
Thrust-zone in the Assynt district of Sutherland, Scotland.’ The fabric of this 
specimen is of particular interest as some of its constituent grains contain 
internally rotated lamellae of the type recorded in experimentally deformed 
dolomite but not hitherto described in naturally deformed rocks. A dynamic 
interpretation of the {0221} twin lamellae and rotated lamellae is made in 
the light of the experimentally established data, and the mechanism of de- 

formation of the rock is tentatively discussed. 


INTRAGRANULAR DEFORMATION OF DOLOMITE 


Inferred from Natural Fabrics.—The first detailed study of dolomite 
orientation in deformed rocks was made by Fairbairn and Hawkes (1941). 
Analyses of the orientation of the lattice and twin lamellae in a number of 
specimens from Montana, Vermont and Ontario were recorded and a possible 
mechanism of deformation was inferred from the data. The authors found 
that {0221} twin lamellae were extensively developed in all the specimens, 
and the orientation of the lamellae in the fabric of the rocks led them to con- 
clude that the sense of gliding in twinning on {0221} was opposite to that for 
twinning on {0112} in associated calcite, They considered that translation- 
gliding on {0001}, which had been produced experimentally by Johnsen 
(1902), might also be of some importance as an orienting mechanism, though 
this would be difficult to prove in view of the difficulty of identifying the 
process by optical means. 

In a comprehensive study of Alpine dolomite-tectonites, Ladurner (1953) 


* The specimen (M14) was collected from a body of crystalline dolomite underlying the 
Moine Thrust at Benmore Lodge near Loch Ailsh. It is one of five similar specimens 
analyzed as part of a detailed structural study of the area. The fabrics of the remaining 
four specimens are described elsewhere (Christie, 1956) and related to the visible mega- 
scopic structures in the rocks of the thrust-zone. 
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recorded the orientation of [0001]-axes in more than 40 specimens from 
various parts of the Alps. On the basis of the preferred orientation of [0001 }- 
axes he classified the rocks into S-tectonites and three types of B-tectonite, 
according to the number of s-planes and symmetry planes which he recognized 
in the fabrics. The similarity between synoptic diagrams showing maxima of 
[0001]-axes in dolomite-tectonites and analogous diagrams for calcite- 
tectonites led Ladurner (1953, p. 290, 296) to assume that the same mechan- 
ism was responsible for producing the preferred orientation in dolomites as 
was postulated by Felkel (1929) for calcite: namely, gliding on {0112}, com- 
bined with rotation of the {0112} planes about the B-axis of the fabric. This 
conclusion is invalidated by the experimental work of Turner, Griggs and 
Heard (1954), who have shown that translation-gliding in calcite occurs, not 
on {0112}, as thought by Felkel, but on {1011}, and at room temperatures 
on {0221}. Moreover, experimental work on dolomite up to the present time 
(see below) gives no indication that gliding on {0112} occurs in dolomite, 
nor is such a mechanism compatible with the geometry of the dolomite lattice. 

Determined by Experimental Deformation of Dolomite.—Early attempts 
to produce twinning experimentally in délomite were unsuccessful but Johnsen 
(1902) recorded translation-gliding in deformed crystals with {0001} as the 
glide plane and one of the a crystal axes as the glide line. More recently two 
glide mechanisms were demonstrated for dolomite (Dover Plains dolomite 
rock) compressed 9.4 percent at 380°C and 3,000 atmospheres confining 
pressure (Turner, Griggs, Heard and Weiss, 1954) : 

l. “twin-gliding on {0221}, the sense of shear being such that the upper 
layers of the crystal lattice are displaced downwards from the upper 
end of the c-axis,” that is, with a negative sense of shear (Turner, 
Griggs and Heard, 1954, p. 897) ; 

2. “translation-gliding on {0001}, with the a-axes as probable glide di- 

rections.” 

Deformation of single crystals of dolomite at room temperature and at 
300°C (Higgs and Handin, 1954) confirms the above conclusions: the crystals 
deformed at 300°C by translation-gliding on {0001} and by twin-gliding on 
{0221} with a negative sense of shear; the latter mechanism was not identified 
in crystals deformed at room temperature. Twin-gliding on {0221} was also 
found by Handin and Fairbairn (1955) in their experiments on Hasmark 
dolomite rock. Bradley, Burst and Graf (1953) have shown, moreover, that 
these two glide mechanisms would account for the progressive changes ob- 
served in the lattice of dolomite after prolonged grinding. 


DYNAMIC INTERPRETATION OF LAMELLAE 

IN NATURALLY DEFORMED DOLOMITE ROCK 
Inferences of two independent kinds, concerning the dynamics of de- 
formation, may be drawn from {0221} lamellae in naturally deformed dolo- 

mite rocks: 

The first is analogous to the dynamic interpretation of twin lamellae in 
calcite marbles, described by Turner (1953). Figure la shows the orientation 
of {0221} twin lamellae in relation to the [0001 ]-axis in dolomite twinned on 
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(0221}. The angle between [0001] and the pole of the lamellae, {0221}, is 
624°. The applied stresses which would be most efiective in causing twin- 
gliding on {0221} with a negative sense are also shown in figure la: a com- 
pression applied in a direction designated C, or a tension in a direction T. 
For either, the coefficient of resolved shear stress on {0221} would have the 
maximum possible value, 0.5. 


Fig. l. a. Projection showing the orientation of [0001] in the host lattice (C,) 
and in the twinned lattice (C,’) in dolomite. C and T are respectively the axes of com- 
pression and tension which give maximum resolved shear stress on {0221} == &.os 
gliding with a negative sense. 

b. Projection showing the orientation of {Q22]} and L9 lamellae in the dolo- 
mite lattice. For fuller explanation, see text. 

c, d. Diagrams illustrating the internal rotation of {Q22]} lamellae, f; and fs, 
to L., and Le, by translation-gliding on (0001) parallel to a:. C., as, a2, as are the crystal 
axes, 


The second kind of inference is drawn from internally rotated {0221} 
lamellae of the type designated L, by Turner, Griggs, Heard and Weiss 
(1954). Following Turner, et al., it is assumed that these have been rotated 
by translation on {0001} parallel to one of the a crystallographic axes. In the 
experimentally deformed Dover Plains rock only one set of L, lamellae were 
found in any grain; but in the specimen here described two sets of L, lamel- 
lae are present in several grains, allowing unique identification of the active 
glide direction. The orientation of these two sets of lamellae in the dolomite 
lattice is shown in figure 1b. Figures lc and d illustrate diagrammatically the 
phenomenon of translation-gliding on {0001} in dolomite. The lattice orienta- 
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tion (C,, a,, a,, ay) is unchanged by gliding. Pre-existing {0221} lamellae, 
f, and f,, are internally rotated about axes which are the intersections of the 
lamellae {0221} with the glide plane {0001}, that is, the a crystal axes. The 
sense of internal rotation about the axes of internal rotation is the same as 
the sense of shear on the glide plane. Following gliding parallel to a, the 
lamellae f, and f, are rotated about a, and a, respectively and so assume the 
new orientations L», and L»,. If three sets of {0221} lamellae are present in a 
grain before the inception of translation-gliding on {0001}, the axes of po- 
tential internal rotation are parallel to the three a-axes of the lattice. But it is 
evident that only two sets of lamellae can undergo rotation by gliding, since 
the axis of potential rotation for the third set is parallel to the glide line. 
Since a, and a, are the axes of internal rotation of L», and L», the only pos- 
sible glide direction is a, and the sense of slip on the glide plane is given by 
the sense of rotation of the L, lamellae. Optimum directions of compression 
and tension which would give maximum resolved shear stress on the {0001} 
glide system so deduced are C, and T,. Such axes of compression and tension 
may be determined for all grains in which two sets of L, lamellae are re- 
corded. For grains in which only one set of L, lamellae is present or is 
accessible for measurement, there is a choice of two possible glide directions 
and there are two alternative orientations of both C and T. For example, if 
only L», lamellae are present (fig. 1b), the axis of internal rotation of L», is 
a, and the active glide line may be either a, or a,; the inferred directions of 
applied stress are C, and T, or C, and T,. 


THE LOCH AILSH DOLOMITE 

Description of the Fabric.—The specimen (M14) described below is a 
pure, massive crystalline dolomite with no trace of foliation or lineation. The 
only impurity in the rock is quartz, in the form of small, isolated grains 
making up considerably less than 1 percent of the rock. The dolomite grains, 
which are more or less uniform in size (average mean grain diameter 
0.54 mms.), show a high degree of post-crystalline strain: all contain well- 
developed {0221} twin lamellae, and a few show undulose extinction. There 
is considerable granulation along grain boundaries, the minute interstitial 
granules comprising approximately 27 percent of the rock (measured by 
means of a point counter in 3 mutually perpendicular sections). There is a 
weak dimensional orientation of the main grains: the majority of the grain 
sections in all three sections are equidimensional, but a limited number in 
each section are slightly elongate parallel to a common direction in the fabric. 

To overcome limitations due to orientation of any single section relative 
to the fabric, the analysis was carried out on three mutually perpendicular 
sections cut from the specimen. Lamellae which are inclined at low angles 
(less than 35°) to a section are not accessible for measurement with the uni- 
versal stage, so that there is a corresponding central “blind spot” in every 
fabric diagram. Since the three {0221} planes in dolomite are inclined to 
each other at approximately 80°, one set of lamellae is inaccessible in most 
grains. It should be noted that L, lamellae can be recorded satisfactorily only 
if the intersection [L,: 0221] can be rotated into parallelism with the micro- 
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scope axis. Consequently, in any section, L, lamellae can be measured only in 
grains with a limited range of orientation. The range of orientation of grains 
in which two sets of L, lamellae can be measured is obviously further re- 
stricted. Thus, in order to obtain a reliable statistical estimate of the orienta- 
tion of L, lamellae in a rock it is necessary to examine at least three sections 
with different orientation. 


Fig. 2. Orientation data for the Loch Ailsh dolomite. All data are plotted on a 
lower hemisphere, equal area projection. Geographic orientation is given by the horizontal 
plane and east (E) and south (S) directions. 

a. 310 [0001]-axes of dolomite. Contours: 1/3, 1, 2, 3, 4, 5% per 1% area. 
b. 300 a@ crystallographic axes in 100 unselected grains. Contours: 1, 2, 3% 
per 1% area. 


c. Poles of 60 twinned (0221) lamellae. Contours: 144, 5% per 1% area. 
d. Poles of {0221} lamellae (points) and associated L, lamellae (arrow- 


heads) in 35 grains. 


Figure 2a shows the preferred orientation of 300 [0001]-axes, 100 
measured in each of three mutually perpendicular sections. A strong similarity 
in the three component partial diagrams (each containing 100 axes) indicates 
a high degree of homogeneity within the field of the specimen, There is a 
single area of concentration, containing two maxima of equal intensity ap- 
proximately 30° apart, and there is some suggestion of a girdle, the axis of 
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which is nearly horizontal in figure 2a. The orientation of the a crystallo- 
graphic axes (figure 2b) is restricted by the strong preferred orientation of 
the [0001 ]-axes: they are disposed in a broad zone normal to the maximum 


PLATE 1 


B 


A. Photomicrogr iph showing the texture of the Loch Ailsh dolomite. Scale line 
represents 1mm. 
B. Single grain showing Ll» lamellae intersected by later {0221} lamellae. The in- 


tersection of the two sets of lamellae is vertical. Scale line represents .lmm. 
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of [0001 ]-axes, and there is a concentration near the periphery of the dia- 
gram (i.e. parallel to the girdle axis of the [0001] diagram). 

As shown in plate 1A {0221} lamellae are strongly developed in all the 
grains. An analysis of the lamellae is given in tables 1 and 2. The majority of 
the grains contain three sets of {0221} lamellae; but it is remarkable that 
there is optically recognizable twinning in only 20 percent of the grains, and 
in these it is generally present on only one set of lamellae. The orientation of 
recognizably twinned lamellae (60 lamellae in 58 grains) is shown in figure 
2c. The poles of the lamellae define a single maximum. 


Tasie 1* 
Percentage of Grains Showing Twinned and Non-twinned {0221} Lamellae 


Number of sets of lamellae per grain Percentage of grains 


3 sets 71% 
2 sets 25% 
1 set 4% 
No sets 0% 


Rotated L, lamellae were recorded in 35 of the grains examined. These 
lamellae (plate 1B) differ from {0221} lamellae in being less sharply defined ; 
they have a discontinuous, granular appearance and their orientation is con- 
sequently more difficult to measure accurately than that of twin lamellae and 
cleavages. The L, lamellae may frequently be recognized by this characteristic 
appearance; but as this is not invariably so, and to avoid errors of identifica- 
tion, they were recorded only where associated with later lamellae having 
rational {0221} orientation (plate 1B). Turner, Griggs, Heard and Weiss 
(1954) record only one set of L9 lamellae per grain in their experimentally 
deformed dolomite rock, but two sets of L, lamellae are present in four of 
the grains analyzed in the Loch Ailsh rock. Figure 2d shows the orientation 
of the L, and associated {0221} lamellae in the specimen; all but two of these 
sets of L, lamellae were measured in the section with the same orientation as 
the diagram. The angle between L, and {0221} varies between 5° and 12°, 


TaBLe 2* 
Percentage of Grains Showing Twinned {0221} Lamellae 
Number of sets of lamellae per grain Percentage of grains 
3 sets 0% 
2 sets 1% 
1 set 
No sets 


The above analyses are based only on grains in which all three planes were accessible 
for measurement with the U-stage. Lamellae parallel to {0221} which are not recogniz- 
able as twin lamellae are described as “non-twinned” to distinguish them from lamellae 


in which the orientation is visibly different from that of the host grain (“twinned lamel- 
lae”), as in Borg & Turner, 1953. 
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with a mean value of 814°, a figure which is comparable with that obtained 
in the dolomite deformed by Turner et al. The sense of rotation of {0221} to 
L, is statistically constant. 

Dynamic Interpretation of the Fabric—The axes of compression and 
tension which would be most effective in causing the observed twinning are 
shown in figure 3a. It is clear that a strong compression parallel to the direc- 
tion C, or a tension parallel to the direction T could account for all the 
observed twinning. Figure 3b shows the glide lines (a) and axes of compres- 
sion (C) and tension (T) deduced for the four grains in which two sets of 
L, lamellae were recorded. For the 31 grains with one set of visible L, lamel- 
lae the 62 possible glide lines are shown contoured in figure 3c (unshaded 


c 


Fig. 3. Dynamic interpretation of data from the Loch Alish dolomite. 

a. Axes of compression (points) and tension (crosses) which would give maxi 
mum resolved shear stress favorable for twinning on the observed twinned {()221]} 
lamellae. 

b. 


Active glide lines (circles) and axes of compression (points) and tension 
(crosses) 


which would give maximum resolved shear stress favorable for translation- 
gliding on {0001} in grains containing two sets of Le lamellae. 
c. 


Synoptic diagram showing maxima of possible glide-lines (a), tension axes 
(T) and « 


ompression axes (C) inferred for 31 grains containing one set of L. lamellae. 

d. Kinematic and dynamic interpretation of the data for the specimen. S, and 
S. are planes defined statistically by {0001}. C,, T; and Ce, T; are stress axes inferred 
from twinned {0221} lamellae and L, lamellae respectively. 


A. im 
x «x *x / T 
7 x/ 
CV /A 
hy 
—— 
A 
d 


Dolomite from the Moine Thrust-zone in North-west Scotland 167 


contours). The possible axes of compression (shaded contours) and tension 
(stippled contours) deduced for the 31 grains are also shown in the figure 
(3c). Maxima for possible glide lines (a) and for axes of compression (C) 
and tension (T) in figure 3c have orientations close to the glide lines and 
axes of compression and tension inferred for grains with two sets of L, 
lamellae (fig. 3b). The axes of compression and tension deduced from twin- 
ning (C, and T,, fig. 3d) are also subparallel to those inferred from the L, 
lamellae (C, and T,, fig. 3d). 

Absence of L, lamellae in a grain does not denote absence of {0001} 
translation, since L, lamellae appear only if {0221} lamellae were present in 
the grain prior to gliding on {0001}. The lattice orientation of grains with 
L. lamellae (figs. da, b) is not obviously different from that of the other 
grains (figs. 2a, b); so that many or all of the grains in the rock may well 
have undergone some degree of translation-gliding on {0001}. Thus the planes 
defined statistically by the orientation of {0001} could represent statistical 
slip-planes in the fabric (S, and S,, fig. 3d). Kinematic a-axes in these slip- 
planes would be given by the preferred orientation of the glide directions in 
individual grains (a). 


Additional orientation data for the Loch Ailsh dolomite. 
[0001 }-axes in 35 grains containing L, lamellae. 
a crystallographic axes in the same 35 grains. Contours: 1, 3, 5% per 1% 


The geometry of the dolomite lattice is such that a compression applied 
subparallel to [0001] will tend to produce twin-gliding simultaneously on two 
or three of the {0221} planes. But if a tensile force is applied normal to 
[0001] in such a way as to give high resolved shear stress on one {0221} 
plane, the resolved shear stress on the other two {0221} planes will be low and 
twin-gliding will take place on only one set of {0221} lamellae. Thus we might 
expect to find two or three sets of twin lamellae in each grain in rocks which 
have been deformed by compression, whereas tension should give a fabric in 
which most grains contain only one set of twin lamellae. The presence of two 
or three sets of {0221} lamellae in the majority of grains in the specimen 
under consideration indicates that these have been produced by compression. 
Although grains containing three sets of {0221} lamellae are very common, 
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only 20 percent of the grains contain optically recognizable twinned lamellae 
and the majority of these have only one such set. This feature, which is rather 
remarkable in view of the considerable post-crystalline deformation of the 
rock, may perhaps indicate a low temperature of deformation, for twinning 
has been produced experimentally only at temperatures above 300°C (Higgs 
and Handin, 1954). 

Although the {0221} lamellae in some grains are closely spaced, there 
is no indication, even where optically recognizable twinning is present, that 
the grains are more than partially twinned. This fact, along with the paucity 
of optically twinned grains, suggests that twinning was of minor importance 
during the final stages of deformation. The evidence of the fabric, then, sug- 
gests that translation-gliding on {0001} was more important than twin-gliding 
on {0221} during the latest stage of the deformation of the rock. 

Orienting Mechanism.—This study yields little evidence as to the me- 
chanism by which the strong preferred orientation of [0001] was achieved. 
It is probable that the intracrystalline structures described and interpreted 
above reflect only the final stage in the deformation of the rock. To produce 
such a strong preferred orientation by gliding of the types discussed would 
require a very high degree of strain, which would result in extreme elonga- 
tion of the grains. Since the grains are predominantly equant, it is likely that 
there was a strong preferred orientation before the visible structures were 
produced. However, the pattern of preferred orientation is so similar to those 
of dolomite tectonites described by other authors (Fairbairn and Hawkes, 
1941; Ladurner, 1953) that there cen be little doubt that it was produced by 
deformation. 

In an aggregate of dolomite crystals with random orientation subjected 
to compression, translation-gliding on {0001} would tend to give a concen- 
tration of [0001 ]-axes close to the axis of compression, whereas twin-gliding 
on {0221} would produce the opposite effect—a migration of [0001 ]-axes 
from the axis of compression. If these are the only two mechanisms of de- 
formation in dolomite, the final pattern of preferred orientation would depend 
solely on which of the two types of gliding was produced with the greatest 
ease. The observed concentration of [0001 ]-axes could have been brought 
about by basal translation with minor twinning on {0221} under compression 
directed parallel to the present [0001] concentration, but definite evidence as 
to the operative mechanism is lacking. 

Calculation of Strain from Rotated Lamellae.—Theoretically it should be 
possible to estimate the post-crystalline strain in a rock from the degree of 
internal rotation of lamellae as a result of gliding according to a known 
system. In view of the extensive post-crystalline strain in the present specimen 
it should be well suited for such consideration (Turner, et al., 1956, p. 1292). 
The strain, ¢, in individual grains is given by the equation: 


So (cota - cotB) 


where a and £ are the angles between the rotated lamellae and the glide- 


siny 
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plane before and after rotation (selected so that B>a); y is the angle be- 
tween the glide line and the axis of internal rotation and So is the coefficient 
of resolved shear stress. For the glide system under consideration y is 60° 
and Sy is given values of 0.4-0.5 (Turner, et al., 1954). For rotations of L, 
ranging from 5° to 12° the calculated strains are 5 to 18 percent. 

There are certain limitations to the acceptance of these figures as a 
measure of the actual strain in the rock. Estimation of strain from the rota- 
tion of L, lamellae assumes that all the deformaiion was achieved by transla- 
tion-gliding on {0001}, and there is evidence in most grains of limited twin- 
gliding. Moreover, many of the grains contain minute dark granules which 
may have resulted from the disruption of rotated lamellae. It is probable, in 
the writer's opinion, that L, lamellae become too diffuse for measurement and 
are eventually disrupted when the strain is much in excess of the values 
recorded above. Thus the strain recorded in the rotation of L, lamellae in 
individual grains is probably less, and perhaps considerably less, than the 


total post-crystalline strain in the rock. 


CONCLUSION 


The Loch Ailsh dolomite has a strongly-oriented tectonite fabric of un- 
known origin, with almost orthorhombic symmetry. The fabric yields a unified 
picture of translation-gliding on {0001} and twin-gliding on {0221}, dating 
from the final stage of deformation of the rock. This late post-crystalline 
deformation probably took place at a comparatively low temperature and 
was produced by a strong compression directed along an axis plunging 50° 
to N315° E. 
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PHASE EQUILIBRIA AT HIGH TEMPERATURES 
IN OXIDE SYSTEMS INVOLVING CHANGES IN 
OXIDATION STATES* 

ARNULF MUAN* 


ABSTRACT. Determination of phase equilibria in oxide systems involving changes in 
oxidation states requires careful control of atmospheric conditions. The desired O2 partial 
pressure may be obtained experimentally by using vacuum techniques or by mixing O2 
with an inert gas in the range 1 to 10° atm. Under more reducing conditions the desired 
O. partial pressure is best obtained by mixing CO. and Hz, in various proportions. The 
resulting O, partial pressure is calculated from known equilibrium data. 

In the description of phase equilibria in such systems, the O, partial pressure of the 
gas phase is an important parameter. Diagrams showing projections of the liquidus sur- 
face with liquidus isotherms and fractionation curves should be supplemented by curves 
showing equilibrium partial pressures of O, of the gas phase. New criteria for derivation 
of paths of equilibrium crystallization must be established. Compositions of condensed 
phases change during crystallization in a manner which is described by straight lines 
pointing towards the O apex of the model used to represent the system. 

Iron oxide containing silicates are the best-known examples of systems where 
changes in oxidation states are involved. The systems Fe—O, FeOQ—Fe.0;—SiO, and 
MgO—FeO—Fe,0,—SiO, are used for the purpose of illustrating principles. 


INTRODUCTION 

\ knowledge of equilibrium as well as kinetic relationships is necessary 
in order to understand and predict chemical reactions. In cases where reaction 
rates are very low, equilibrium data are of relatively little use for an evalua- 
tion of industrial processes, while they still may be of great importance for 
an understanding of geological processes, where the time factor is less critical. 
As reaction rates increase, equilibrium data become increasingly important. 
Although few, if any, industrial processes operate under perfect equilibrium 
conditions, equilibrium data have their great value in defining limits for the 
reactions taking place. 

It is customary, for classification purposes, to distinguish between 
homogeneous and heterogeneous equilibria. Actually, a single homogeneous 
phase can not be isolated; it is not possible in practice to contain or support 
it without introducing one or more other phases. However, where the main 
reactions take place within the homogeneous part of the system, such as for 
instance in a gas phase, the equilibrium may be considered homogeneous. 
Heterogeneous equilibria, on the other hand, are those in which mass transfer 
across phase boundaries is an important part of the process. Such equilibria 
invariably involve condensed' phases, either liquid or solid or both. The con- 
ditions governing the coexistence of phases at equilibrium can be derived 
from classical thermodynamics, and are commonly expressed in the law re- 
ferred to as Gibbs’ Phase Rule. 


+ c+ 2 


Here P designates the number of phases existing together at equilibrium, F 
is the number of degrees of freedom and C the number of components. The 
dependence of phase relationships on the variables of the system is best il- 


* Contribution No. 57-4 from College of Mineral Industries, The Pennsylvania State Uni- 
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lustrated in the geometrical constructions known as phase diagrams. Among 
the variables influencing the phase assemblage may be temperature, pressure, 
chenical composition, electrical field, magnetic field etc. It has been found 
expe~imentally that the latter two may usually be ignored. Furthermore, in 
most cases reactions are carried out at a total pressure of 1 atm.’, leaving 
temperature and chemical composition as the main variables. The dependence 
of phase relationships on these variables is what is usually illustrated in a 
phase equilibrium diagram. 

A very substantial amount of phase equilibrium work on oxide systems 
has been described in the literature during the last decades.’ In reviewing 
these studies one is impressed by the large amount of data available on 
systems made up from oxides of noble gas type ions, such as for instance 
SiO,, Al,O,, MgO, CaO. The great interest in such oxides is, of course, well 
founded. These oxides are among the most abundant constituents of the 
earth’s crust, and a knowledge of systems made up from these components is 
important in geological sciences as well as in technology. There are, however, 
a large number of other oxides, as well as other compounds, which are equally 
important, and becoming more so as new materials are developed, Many of 
these materials involve components made up from ions of transition elements 
where changes in oxidation states not only are possible but indeed very com- 
mon. The study of phase equilibria in such systems is somewhat more com- 
plicated than the study of the systems just mentioned, because the gas phase 
plays an important role in the equilibria. The present paper discusses methods 
for studying and describing such systems. 


METHODS FOR DETERMINING PHASE EQUILIBRIA IN OXIDE SYSTEMS 

Ordinary Dry Silicate Systems.—In working out and describing phase 
equilibria in “ordinary dry” silicate systems‘ it is customary to ignore the 
vapor phase and treat the equilibria in terms of the condensed phase rule 

+? C + 1 

where P indicates the number of condensed phases present, while the other 
symbols have the same meaning as explained before. No serious error is in- 
troduced by doing this as long as the vapor pressure of each constitutent is 
low. Such systems are usually studied in an open quench furnace in air 
atmosphere. As an illustration, consider the system CaOQ—MgOQ—SiO, at 
liquidus temperatures. At equilibrium there is a certain partial pressure of 
each constituent in the gas phase coexisting with the condensed phases. These 
partial pressures vary with temperature and composition of the liquid, each 


assemblage of condensed phases requiring, at equilibrium, a very definite 


composition of the gas phase. During the time required for equilibrium to be 
* The term condensed phases is used here and throughout the remainder of the paper to 
designate liquid and solid phases as distinguished from a gas phase. 

* Pressure is, of course, an extremely important variable in geological sciences, as evi- 
denced by the large amount of research work presently going on in the field of hydro- 
thermal synthesis at high pressures and elevated temperatures. 

* See for instance: Levin, E. M., McMurdie, H. R. and Hall, F. P., 1956, Phase diagrams 
for ceramists: Columbus, Ohio, Amer. Ceramic Soc. 
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attained among condensed phases the equilibrium composition of the gas 
phase over the charge may or may not be obtained, depending to a large 
extent on the degree or circulation of air through the furnace. Whether or not 
this true equilibrium composition of the gas phase over the sample is reached, 
is, however, not important in practice. This is because the vapor pressure of 
each constituent is so low® that the sample may build up, or try to build up, 
the equilibrium composition of the gas phase without changing appreciably 
the total composition of the condensed phases. Equilibria among the latter 
can therefore be studied without controlling the atmosphere, and the phase 
relationships can be discussed in terms of the condensed phase rule. 

Systems Involving Changes in Oxidation States.—In cases where elements 
of changing oxidation states are involved in the equilibria, the situation is 
usually more complicated. Consider an arbitrary element Me forming two 
oxides MeO, and MeO,. Here the element Me occurs in two different states of 
oxidation, Me** and Me*’, respectively. Suppose that these oxides are dissolved 
in a condensed phase in equilibrium with a gas phase, The following reaction 
is of critical importance: 


MeO, soln, (MeO, ) in soln 


It will be seen that a change in O, pressure causes a change in the relative 
proportions Me**, Me®’ in the condensed phase, and vice versa. Hence in treat- 
ing systems involving such components the partial pressure of O, must be 
taken into consideration in addition to the vapor pressures of the components 


themselves. As long as all these pressures are very low (of the same order of 
magnitude as for CaO, MgO and SiO,.), such systems can be studied and 
described in a manner similar* to the treatment of CaQ—MgO—SiO,. The 
vapor pressures of the oxides themselves as a function of temperature are 
predetermined. For a large number of important oxides, such as for instance 
iron oxides, these pressures are very low up to liquidus temperatures, The 
equilibrium O, pressure, however, is a function of the ratio MeO, MeO, in 
the condensed phases, increasing rapidly as this ratio increases. If the MeO,/ 
MeO, ratio of the condensed phases in a mixture is such that the equilibrium 
QO, pressure of the gase phase is low, phase equilibria for this mixture can be 
studied in a vacuum system or in a purified inert gas atmosphere. The basis 
for the use of this method is the slow exchange of O, between the gas phase 
and the condensed phases. As an illustration, assume that the MeO,/MeO, 
ratio in the liquid phase at the temperature considered is such that the cor- 
responding partial pressure of O, is 10~° atm. Assume furthermore that a 
completely “perfect” inert gas is used, with no O, present in the gas initially. 
The mixture when heated in a furnace through which this inert gas flows will 
try to build up the equilibrium O, pressure of 10~-* atm. The amount of 
* The term “ordinary” is used in this paper to designate oxide systems involving only ions 
of noble gas type, as distinguished from systems containing transition elements. The term 


“dry” is used to indicate that no water is present, as distinguished from hydrothermal 
systems. 


* The vapor pressures of CaO, MgO and SiO, at 1400°C in air are approximately (Brewer, 
1953): CaO, 10°” atm.; MgO, 10° atm.; SiO., 10° atm. 
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oxygen required to build up this low pressure in a moderate volume of gas 
(the inert gas flows through the furance at a moderate speed) is very small 
9 

( —e . 32 1.4 X 10-° g oxygen per liter of gas). Hence the MeO, 

MeO, ratio of the condensed phases is changed very little, and for practical 
purposes the composition of the condensed phases may be considered con- 
stant, If the inert gas is not “perfect”, the situation is similar. Assume for the 
sake of argument that the inert gas (or the vacuum system) contains a small 
amount of oxygen corresponding to an O, pressure of 10~' atm. In order to 
establish equilibrium between the gas phase and the condensed phases the 
latter phases must pick up O, from the gas phase. However, the difference 
10-* — 10~° is very small, and again an enormous volume of gas would be 
required to change significantly the relative amounts MnO,/ MnO, and hence 
the total composition of the condensed phases (~1.4 X 10~* g oxygen per 
liter of gas). For all practical purposes, therefore, the total composition of 
the condensed phases may be considered constant. This is the basis for the 
method used by Schuhmann, Powell and Michal (1953) in their study of a 
part of the system FeOQ—Fe,0,—SiO.. They made up mixtures in the low 
FeO part of the system by weighing up the three components FeO, FeO 
and SiO, in required proportions. The mixtures were heated in platinum con- 
tainers in an inert gas atmosphere, and the results interpreted on the basis of 
the assumption that the total composition of the condensed phases did not 
change during the period of time necessary to attain equilibrium among the 
condensed phases. 

An alternative method for studying such systems under reducing condi- 
tions is the one used by Bowen and Schairer (1932) in their investigation of 
phase equilibria in iron silicates. Their work, which was started 25 years ago, 
was the pioneering effort in the field and represented the first real break- 
through in experimental studies of iron oxide containing systems. The method 
consists in examining equilibria among condensed phases in mixtures con- 
tained in Fe crucibles in a purified inert gas atmosphere, such as for instance 
N.. The partial pressure of O, of the gas phase under these conditions is that 
corresponding to the equilibria 


l 
Fe.O ), n melt) 7 9 0, 


(g) 


and Fe.03;in meit) + 3FeO; in meit) 
at any particular temperature. In the binary system Fe — O only one degree 
of freedom exists in the presence of the phases above (liquid, metallic Fe, 
gas), and the O, pressure is a function of temperature alone. In a ternary 
system, such as for instance FeQ—Fe,0,—SiO,, one additional degree of 
freedom exists with the same phase assemblage, and the equilibrium O, pres- 
sure is a function of temperature as well as the amount of one of the com- 


* Air atmosphere can not be used for the experimental study of such equilibria. This point 
will be discussed in more detail below. 

* These authors were fully aware of the composition changes caused by iron going into 
the platinum of the container, according to the reaction 3FeO Fe + FesOs. However, 
in most instances this reaction was not extensive enough to change the results significantly. 
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ponents of the liquid phase. If wiistite (FeO)* is one of the phases present, 
this O, pressure is equal to the equilibrium pressure at which FeO and Fe 
coexist in equilibrium in the binary system Fe — O, at the temperature in 
question. (For data on this equilibrium, see papers by Darken and Gurry 
(1945, 1946). If wiistite is not present as a phase, at temperatures below its 
melting point (1371°C), this O, pressure is lower than the pressure cor- 
responding to the FeO /Fe equilibrium. 

At considerably higher levels of O, pressures, equilibrium may be 
reached between an O, containing gas phase and the condensed phases with- 
in a reasonable period of time. As an example, one liter of O, gas of 1 atm. 
pressure corresponds to 1.43 g oxygen. Hence a moderate rate of flow of gas 
through the furnace in this case will supply enough oxygen to change the 
MeO, MeO, ratio of the condensed phases to the equilibrium value within a 
reasonable period of time. It is therefore possible in practice to study equili- 
bria in oxide systems, involving a gas phase at 1 atm. O,, and in air (cor- 
responding to po, 0.21 atm.). If the O, pressure of the gas is decreased, 
either by using a vacuum or by mixing with an inert gas, the volume of gas 
necessary to bring about the change in composition of condensed phases re- 


quired for equilibrium increases rapidly (proportional to — -) . An O, 


pressure of approximately 10~* atm. is a practical lower limit for studying 
the equilibria when oxygen is present only as O, molecules in the gas phase. 

It is possible, however, to obtain equilibrium between the gas phase and 
the condensed phases also below these levels of O, pressures by using indirect 
methods of creating the desired O, partial pressures. The most convenient 
way of doing this is to use an oxygen-containing gas which upon heating in 
the furnace decomposes to another gas species under liberation of O,. Gases 
suitable for this purpose are H,O and CO,, which decompose upon heating 
according to the equations: 


2H.0 2H. + O., and 
2CO, 2CO + O.,, respectively. 


The equilibrium constant for each of these reactions as a function of tem- 
perature is well known. If CO, or H,O is used alone, and the total pressure 
is kept constant, the O, pressure is a function of temperature alone. This 
means that O, pressure and temperature can not be varied independently if 
the pure gases are used, and this of course is a serious drawback. If, however, 
some H, is mixed with the H.O or some CO with the CO.,® the dissociations 
of H.O or CO, are suppressed to a degree which can be controlled by varying 
the mixing ratios of the gases. The O, pressure now is a function of 
temperature as well as H,O/H, or CO,/CO mixing ratios. O, pressure and 
* Wiistite is a defect structure with a considerable variation in the Fe/O ratio (Darken 
and Gurry, 1945, 1946). Although the stoichiometric composition FeO can not be realized 
in practice, the formula for wiistite is usually written as FeO for the sake of convenience. 
* These gases are interchangeable; CO may be mixed with HO and Hz with CO.. The 
interdependence of the equilibrium ratios HsO/Hs, CO./CO, CO./Hs and H:O/CO at a 


certain O, pressure is expressed by the equilibrium constant for the “water gas equation” 
CO, + Hs + H,O. 
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temperature can therefore be varied independently. The advantage of using 
mixed gas atmospheres rather than vacuum or inert gas is the presence of 
oxygen not only as Q, but in bound form (H,O or CO.) as well, At tempera- 
tures high enough to make the gas reactions rapid, a mechanism is available 
by which the O, exchange between the gas phase and the condensed phases 
can take place within a reasonable period of time in spite of the low partial 
pressure of O,. As soon as O, is removed from the gas phase by reaction with 
the condensed phases, the equilibrium amount is restored by the decomposi- 
tions 2H.O 2H.O + O, or 2CO, 2CO + O., whatever the case may 
be. It is emphasized that such methods for equilibrium studies can be used 
only at temperatures sufficiently high to make the gas reactions rapid. 
REPRESENTATIONS OF PHASE RELATIONSHIPS 

The System Fe — O.—Systems containing iron oxides as components have 
been studied more thoroughly than any other oxide systems in which changes 
in oxidation states are involved. These will be used as examples in the fol- 


lowing. 
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Fig. 1. Diagram showing stability relationships among the phases hematite (Fe.O;), 
magnetite (FeO.Fe.0,), wiistite (FeO), metallic iron (Fe), liquid oxide and liquid iron 
as a function of temperature (in °C) and O, pressure (in atm.). Heavy solid lines are 
boundary curves separating the various phase areas labelled in the diagram. Light dash- 
double dot lines are lines of equal px,o/px, ratios of the gas phase at 1 atm. total pres- 
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A thorough understanding of phase equilibria in the system Fe ~ O is 
essential as a background for discussing more complicated multicomponent 
iron silicates. The most recent and complete reports on the system are the 
excellent papers in 1945 and 1946 by Darken and Gurry (1945, 1946). These 
authors summarize results of earlier work in addition to providing a large 
amount of new data. The following discussion is based on data contained in 
their papers as well as data available in a recent compilation of thermody- 
namic constants by Coughlin (1954). 

The importance of controlling the O, pressure was emphasized in the 
last section. The diagram shown in figure 1 is intended to bring out this 
point more clearly and in a quantitative way for the Fe—O system. In this 
plot phase equilibrium relationships among metallic iron (Fe), wiistite 
(FeO), magnetite (FeO.Fe,0,), hematite (Fe,O,), liquid iron and liquid 
oxide are illustrated as a function of temperature (°C) and partial pressure 
of O,. The diagram has been constructed from data explained in detail later 
in connection with figure 2. The following rules should be observed when 
using the diagram in figure 1. In an area where one condensed phase exists 
in equilibrium with a gas, the system possesses 2 degrees of freedom. Hence, 
in a field such as for instance that of wiistite both temperature as well as O, 
pressure may be changed without causing any phase changes. Along a bound- 
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Fig. 2. Diagram presented in order to illustrate stability relationships among iron 
oxides and metallic iron as a function of temperature (in °C) and ratios of H,O/Hz: of 
the gas phase in equilibrium with the condensed phases. Heavy solid lines are boundary 
curves separating the various phase areas labelled in the diagram. Light dash-dot lines 
are Q, isobars, The heavy dash-cross line shows equilibrium ratios pu.o/pa, of pure 


water vapor as a function of temperature at a total pressure of 1 atm. 
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ary curve, however, where 2 condensed phases coexist in equilibrium with a 
gas phase, the system has only | degree of freedom. Hence, if temperature is 
chosen, the O, pressure is fixed and vice versa. It will be noticed that the O 
pressures of the gas phase in equilibrium with the various condensed phases 
in the system vary over a very large range, as temperatures vary between 
room temperature and liquidus temperatures. This of course is directly related 
to the large heat effects associated with each of the reactions represented by 
the boundary curves. ( = an , Where AH” is the heat of re- 
dT RT? 

action per mole of O, for each reaction, such as for instance: 6Fe,0 

iFe,0, + O,). Oxygen pressures at the lowest temperatures shown in 
figure 1 are so low as to have little or no physical meaning. Physically more 
meaningful are ratios px.o/px, or peo,/peo of the gas phase in equilibrium 
with the various oxides. Lines representing equilibrium values of px,.o/px, 
ratios are plotted as dash double dot lines in figure 1. It is seen that these 
lines fall close to and are nearly parallel to the boundary curves separating 
the various phase areas. More “spread out” and readable diagrams are ob- 
tained by plotting px,o/px, and peo,/peo ratios of the gas phase in equilib- 
rium with the various condensed phases as a function of temperature, as 
shown in figures 2 and 3, respectively. The diagram in figure 2 has been con- 
tructed from the following data: The boundary curve along which metallic 
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Fig. 3. Diagram illustrating stability relationships among iron oxides and metallic 
iron as a function of temperature and ratios of CO./CO of the gas phase in equilibrium 
with the condensed phases. Heavy solid lines and light dash-dot lines have the same mean- 
ings as lif ast lagran The heavy dash-double cross line indicates ratios Peo./Pco in 
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iron and wiistite or liquid oxide coexist in equilibrium with gas (curve II), 
as well as part of the curve where magnetite, wiistite and gas exist together in 
equilibrium (curve III), have been plotted from data obtained by direct 
measurements of H,O/H, ratios of the gas in equilibrium with these phases, 
as reported in the literature. The lower temperature part of curve II (600- 
1000°C) was constructed from the data of Emmetts and Shultz (1930), 
while the upper region (1200-1515°C) was plotted from the data of Chipman 
and Marshall (1940). The low temperature part of curve III (600-800°C) 
was drawn on the basis of the data of Emmetts and Shultz (1930). The upper 
part of this curve (1096-1388°C) was constructed on the basis of data re- 
ported by Darken and Gurry (1946), These authors used CO,/CO and CO,/ 
H. gas mixtures rather than H.O/H, to study the equilibrium described by 
this line. Their data have been converted to values of H,O/H, ratios by use 
of the equilibrium data for the reactions 


2H.0 = 2H, + 0, 


as tabulated in a recent publication by Coughlin (1954). The liquidus curve 
V, as well as that part of the hematite-magnetite boundary curve IV falling be- 
tween 1100 and 1450°C, likewise, have been constructed from the data of 
Darken and Gurry (1946) by conversion to H,O/H, ratios as described 
above. No experimental data are available for the boundary curve magnetite- 
metallic iron (curve I) and the hematite-magnetite boundary curve below 
1100°C (curve IV). These parts of the curves have been constructed on the 
basis of calculated values obtained from the standard free energy data tabu- 
lated by Coughlin (1954). For checking purpose, curves were also constructed 
from calculated values in regions where direct measurements have been made 
(curves II, III, part of curve IV). The calculated curves are in good agree- 
ment with the experimental curves except for the magnetite-wiistite boundary 
(eurve III). The discrepancy in the latter case is not surprising. The data 
listed in the tables of Coughlin refer to wiistite of composition Fe, .,-O, while 
the wiistite in equilibrium with magnetite varies in composition in the range 
Feo 9500 to Feo cso0 (see phase diagram of Darken and Gurry (1945, 1946) ). 
Light dashed lines show partial pressures of O,. The heavy dash-cross line in 
the upper part of the diagram gives equilibrium values of px.o/px, of pure 
water as a function of temperature. The position of this curve relative to the 
boundary curves shows the stability ranges of the various iron oxides in pure 
water vapor. It will be noticed that at temperatures below approximately 
1070°C hematite is the stable phase in equilibrium with water vapor. At tem- 
peratures between 1070°C and 1595° magnetite is the equilibrium phase in a 
water vapor atmosphere, and above 1595°C a liquid oxide phase is stable 
under these conditions. It should be pointed out that neither wiistite nor 
metallic iron are equilibrium phases in water vapor at any temperature at a 
total pressure of 1 atm. In order to study equilibria in regions of the dia- 
gram falling below this dissociation curve for pure H,O it is necessary to mix 
H. (or some other reducing agent) with the H,O in amounts that may be cal- 
culated from the diagram in figure 2. One important relationship should be 


2CO, = 2CO + O, 
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made clear at this point: The values of H,O/H, ratios plotted in figure 2 are 


the equilibrium ratios at any particular temperature. In order to study ex- 
perimentally equilibria among the iron oxide phases it is usually necessary to 
mix the two gases in fixed proportions at room temperature and heat this gas 
mixture to whatever temperature is desired for the study of the equilibria 
among the condensed phases. These two ratios (original mixing ratio at 
room temperature and equilibrium ratio at the elevated temperature) are 
not the same. They are related to each other in the following way: Let pu,o 
and pw.«> be the partial pressures of H,O and H,, respectively, in the starting 
gas mixture at room temperature. Upon heating, H.O decomposes according 
to the equation 
2H.0 = 2H, + O, 
At equilibrium we have the relationship 
IN 
p 
where p etc. are the partial pressures of the gases at equilibrium at the 
high temperature, and K is the equilibrium constant at the same temperature. 
Let us call p x. Then we have the following relationships: 


| ) 


and likewise 


These approximations are valid in all regions of the diagram in figure 2 ex- 


cept in close proximity of the line for dissociation of pure water, For pure 


water the following relationships hold: 
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At moderate or low temperatures, where K is small, x< <1, and the following 
approximation is valid: 


Hen e 


In order to study equilibria above the dissociation curve for pure water 
(fig. 2), on the other hand, it is necessary to mix O, with the H,O. At mod- 
erate or low temperatures, where the dissociation of water is small, this 
amounts to practically the same as mixing O, with an inert gas; the partial 
pressure of O, in the mixture is rcactically the same as the proportion of O, 
mixed with the gas. 

It will be observed that some of the boundary curves in figure 2 have a 
positive slope, some have a negative slope, and one (the Fe—wiistite and Fe— 
liquid oxide boundary) is practically vertical. The latter case illustrates a 
situation in which the ratio px.o/ps, of the gas in equilibrium with the two 
condensed phases remains essentially constant over a large temperature in- 
terval (1000-1600°C). The reason for this is that the heat of reduction of iron 
oxide to Fe is balanced by an almost identical heat effect of opposite sign, 
the oxidation of H. to H.O: 


FeO + H, Fe + H,O 


Therefore, AH” for this reaction is very small, and the equilibrium constant 
(ratio H,O/H,) is almost independent of temperature. Curves with positive 
slope (for example the Fe—Fe,O, boundary curve) result when the heat of 
reduction of the oxide to the next lower oxide or to metallic iron is larger 
than the heat of oxidation of H, taking part in the reaction. Conversely, the 
negative slope of the hematite-magnetite boundary curve illustrates an example 
of an overall exothermic reaction. 

The diagram in figure 3 shows equilibrium ratios poo,/pes of the gas 
phase in stable equilibrium with condensed phases in the system Fe—O, as a 
function of temperature. The diagram has been constructed from the same 
data as explained in connection with figure 2. It will be observed that the two 
diagrams in figures 2 and 3 are almost identical. The reason for this is that 
the dissociation constants of H,O and CO, are closely similar, the former 
being slightly higher at temperatures above 800°C, the latter being slightly 


P, x 
Pu 2x 
x* (2x)? 
K \° 
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higher at lower temperatures. The various lines have the same meanings as in 
the last diagram (fig. 2) with pco./pco ratios replacing the p=,0/ps, ratios. 
An additional line of great importance in the present case is the dash-double 
cross line illustrating peo./pco ratios in equilibrium with graphite according 
to the equation 


= Cre) + CO, 


2 (g) 

This line indicates limiting ratios of peo./peo at various temperatures. If, 
at a chosen temperature, CO, CO ratios of the gas phase are decreased from 
a value represented by a point in the upper part of the diagram, graphite 
starts depositing, with CO, CO ratio remaining constant, when the latter ratio 
reaches a value corresponding to the point on the dash-cross curve for this 
temperature. The reaction 2CO,, C,.) + CO, .., is strongly exothermic, 
as evidenced by the steep negative slope of the dash-cross curve in figure 3. 
In order to obtain a gas with low CO,CO ratios, it is therefore necessary to 
go to high temperatures. It will be noticed that the dash-double cross line 
intersects the Fe,O,—FeO and the FeOQ—Fe boundary curves at approximate- 
ly 660 and 720°C, respectively. These are the lowest temperatures at which 
wiistite and metallic iron, respectively, can exist in equilibrium with C - O 
gases. This relationship of course is of fundamental importance for an under- 
standing of the blast furnace process, and has some important implications in 
certain fields of geochemistry. During the last couple of years geochemists 
have taken a great deal of interest in equilibrium relationships of carbonates 
at high CO, pressures. (See for instance papers by Harker and Tuttle, 1956; 
Harker and Hutta, 1956). One of the carbonates of great interest in geo- 
chemistry is siderite (FeCO,). The graph in figure 3, although constructed 
for a total pressure of | atm., may be used to judge the possibility of studying 
experimentally the reaction 


FeCO,,, = FeO,-) + CO; 


(gs) 
It is evident from the diagram that a considerable proportion, the exact 
amount depending on temperature, of CO must be mixed with the CO, in 
order to prevent oxidation of FeO to Fe,O,. (See curve III in fig. 3), The 
disproportionation of CO to CO, and C, represented by the dash-double cross 
curve, however, puts a limit on the CO/CO, ratios that can be reached. The 
intersection between this curve and curve III at 660°C defines a limiting tem- 
perature below which FeO can not exist as an equilibrium phase in the 
presence of C and O at a total pressure of 1 atm. Higher pressures will cause 
a shift in the dash-double cross curve towards the right, while the boundary 
curve III may be expected to be comparatively little effected by an increase 
in pressure. The lowest temperature at which siderite may be decomposed to 
give wiistite as an equilibrium phase, therefore, is higher than 660°C at high- 
er pressures. The metastable equilibrium FeCO, = FeO + CO, may be 
calculated, below these temperatures, from experimental equilibrium data on 
the reaction 3FeCO Fe,0, + CO + 2C0,, if such become available, 
combined with thermodynamic data referred to above (Darken and Gurry, 
1945, 1946; Coughlin, 1954). 
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Fig. 4. Diagram to show phase equilibrium relationships in the system Fe—O, 
mainly after Darken and Gurry (1945, 1946). Heavy solid lines are boundary curves 
separating phase areas labelled in the diagram and light dash-dot lines are O. isobars 
(in atm.). Dashed curves in upper right part of the system are inferred boundary curves 
in a region where experimental data are lacking (po, >1 atm.). 


The diagrams presented in figures 2 and 3 have one serious drawback: 
They do not show compositions of condensed phases. It is of course possible 
to draw a phase diagram for the system Fe—O showing compositions of con- 
densed phases. This indeed has been done; the most complete diagram is that 
presented by Darken and Gurry (1945). They plot the diagram in the usual 
way with temperatures (in °C) along the vertical axis and composition (ex- 
pressed in weight © O) along the horizontal axis. Their diagram also has 
one serious drawback, however: It does not show compositions of the gas 
phase in equilibrium with the condensed phases. The diagram presented in 
figure 4 of the present paper is intended to illustrate phase equilibria in the 
system Fe—O in a manner which eliminates the drawbacks mentioned above. 
The diagram is essentially a reproduction from Darken and Gurry (1946), 
the axes being temperature and composition. In addition to the boundary 
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curves shown in their diagram a family of curves representing equal values 
of partial pressures of O, of the gas phase in equilibrium with the condensed 
phases is shown in figure 4. It would be possible, similarly, to construct a 
family of curves representing equal px,o,//px, or pco,/peo ratios of the equilib- 
rium gas phase. 

The relationships among the various types of representations should be 
kept clearly in mind in the use of the diagrams in figures 2-4. In figures 2 
and 3 each area corresponds to a situation where one condensed phase, as 
labelled in the diagram, exists in equilibrium with a gas. Boundary curves in 
the same two diagrams represent situations where two condensed phases and 
a gas coexist in equilibrium, In figure 4, on the other hand, areas in the dia- 
gram may represent situations where either one or two condensed phases 
coexist in equilibrium with a gas phase. These observations are not conflict- 
ing. In figures 2 and 3 there is no axis to represent compositions of the con- 
densed phases. The two condensed phases of different compositions existing 
in equilibrium along a boundary curve are therefore recorded in these dia- 
grams by one point, namely that point representing the combination of 
temperature and H,O/H, (or CO,/CO) ratio at which these phases exist in 
equilibrium. In fig. 4, however, the axes have been chosen so that compositions 
of condensed phases are shown. Because the compositions of the condensed 
phases in equilibrium are different, the two phases are represented by sep- 
arate points in this diagram. It is as if the boundary curves in figures 2 and 
> were “pulled out” to form areas in figure 4. In figures 2 and 3 adjacent 
phase areas are separated by boundary curves, in figure 4, correspondingly, 
each one-phase area (one condensed phase) must be limited on either side by 
an area where two condensed phases coexist in equilibrium with a gas phase. 
Because the system possesses two degrees of freedom when one condensed 
phase exists in equilibrium with a gas, both temperature as well as O, pres- 
sure or H.O/H, (or CO./CO) ratio may vary in one-phase areas in figures 
2-4. Therefore the lines representing O, isobars in figure 4 run diagonally 
across such areas. In a two-phase area (two condensed phases), on the other 
hand, only one degree of freedom exists and po, as well as H,O/H, (or 
CO.,/CO) ratios in the gas phase are fixed if the temperature is chosen. The 
O, isobaric lines therefore must be horizontal lines parallel to the composition 
axis in a two-phase region. 

Ternary and Quaternary Iron Silicates——The methods used for describ- 
ing phase equilibrium relationships in the binary system Fe—O can be ex- 
tended into the more complicated ternary and quaternary systems. As an 
illustration of the former, consider the system FeQ—Fe,0,—SiO,, which is 
part of the ternary system Fe—Si—O. In the following we will focus our at- 
tention on the liquidus temperature region, In figure 5 is shown the usual 
plot of the projected liquidus surface with boundary curves and liquidus 
isotherms (Muan, 1955). When one crystalline and one liquid phase coexist 
with a gas phase, the system possesses two degrees of freedom. If temperature 
and liquid composition are chosen, O, pressures are fixed. It is therefore 


possible to draw a set of curves representing equal partial pressures of O, 
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Fig. 5. Phase equilibrium diagram at liquidus temperatures for the system FeQ— 
Fe,O,—SiOs, after Muan (1955). Light lines are liquidus isotherms at 100°C intervals, 
and heavy lines are boundary curves, with arrows pointing in directions of falling tem- 
peratures, The tridymite-cristobalite boundary curve is dashed. Medium lines with stip- 
pling on one side indicate limit of two liquid region. 


(O, isobars) of the gas phase in equilibrium with condensed phases along the 
liquidus surface. These lines are extensions into the triangle representing the 
FeO—Fe,0,—SiO, system of points of intersection between O, isobars and 
the liquidus curve in the Fe—O system (fig. 4). An example of such lines is 
illustrated in figure 6 (Muan, 1955). Similarly, one may draw lines repre- 
senting equal pu.o/px, (or peo,/peo) ratios for points along the liquidus 
surface. Figure 7 gives an illustration of this for the system FeOQ—Fe,0,— 
SiO, (Muan, 1955).'° It will be noticed that the lines of equal peo. px, ratios 
are nearly parallel to the boundary curve between the fields of wiistite and 
metallic iron. This is in accordance with figure 3, where it was shown that 
this boundary curve is almost vertical (corresponding to constant peo., peo 
ratio)*” over a large temperature interval. 

’ This diagram shows lines representing equal pco./pu., ratios rather than peo./co 
ratios. However, the numerical values of these ratios are almost identical because the 
equilibrium constant for the “water gas reaction” CO, + H, CO + H.O is not sig- 
nificantly different from 1 in the temperature interval considered. 


185 
20, 
30, ro 
40 CRISTOBALITE so 
80, 
40 
f 
° ~—S 
g 
os ~ ~ 
wf =~ 
MAGNETITE 
| aoe 
FeO 20 » 40 so 


{rnulf Muan—Phase Equilibria at High Temperatures in 


$0 60 FeO Fe,0, 


WEIGHT PER CENT 


Fig. 6. Phase equilibrium diagram for the system FeO—Fe.0;—SiO,. Heavy lines 
are boundary curves and dash-dot lines are lines of equal O, pressures for points on 


liquidus surface. Lines with stippling on one side indicate limit of two liquid region. 
(Diagram after Muan (1955) ). 


Liquidus temperatures were chosen for illustration in the present case. 
One may of course construct O, isobars or lines representing equal px,o/px, 
(or peo. peo) ratios at other temperatures as well, above or below liquidus 
temperatures. 

In order to illustrate the relationships of compositions of crystalline 
phases in equilibrium with liquid and gas, the above diagrams are supple- 
mented by the fractionation curves shown in figure 8. Fractionation curves 
are constructed from conjugation lines in such a way that the conjugation 
lines are tangent to the fractionation curves. (A conjugation line joins the 
composition of a liquid with the composition of the crystalline phase with 
which it is in equilibrium). The composition of a crystalline phase in equilib- 
rium with a liquid, therefore, may be derived from the diagram of figure 8 as 
the intersection with the join FeOQ—Fe,O, of the tangent to the fractionation 
curve at the point representing liquid composition. 

In the case of a quaternary system the principles are the same, but it is 
much more difficult to illustrate the relationships geometrically, Take as an 
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7. Phase equilibrium diagram for the system FeO—Fe.0,—SiO:. Heavy lines 
are boundary curves and dash-double dot lines are lines of equal pco./'pu, ratios of the 
gas phase in equilibrium with liquids at liquidus temperatures. Lines with stippling on 


one side indicate limit of two liquid region. (Diagram after Muan (1955) ). 


example the quaternary"' system MgO—FeO—Fe,0,—SiO, (Muan and 


Osborn, 1956). For our purpose it is most convenient to consider this system 
as formed from FeOQ—Fe,0,—SiO, by the addition of MgO. If, as before, we 
focus our attention on the liquidus temperature region, the O, isobars or the 
lines representing equal pco,/ps, ratios along the liquidus surface in the 
ternary system FeQ—Fe,0,—SiO, are starting lines for irregularly curved 


surfaces (“oxygen isobaric surfaces”) extending through the tetrahedron 
representing the Mg0—FeO—Fe,0,—SiO, system. These surfaces are difficult 
to illustrate in perspective drawings of the tetrahedron, An indication of their 
location is obtained in the following way: Phase equilibria at an O, pressure 
represented by such a surface are illustrated by projecting points in the sur- 
face into a chosen plane cutting through the tetrahedron, preferably a plane 
located close to the surface in question. The actual location of composition 
points of the mixtures at liquidus temperatures, and hence the location of the 
“ The system MgO—FeO—Fe,0,—Si0, is part of the quaternary system Mg—Fe—Si—O 
but is not in itself a completely quaternary system because metallic iron appears as a 


phase in certain composition-temperature ranges. For our purposes this may be disre- 
garded, and the system will be considered quaternary. 
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Fig. 8. Diagram of the system FeO—Fe.O,.—SiO, presented in order to illustrate 
relationship of compositions of liquid and crystalline phases in equilibrium, after Muan 
(1955). Heavy lines are boundary curves and medium dash lines are fractionation curves. 
Lines with stippling on one side indicate limit of two liquid region. 


isobaric surface at liquidus temperatures, is indicated by lines in the pro- 
jected diagram representing equal Fe,O,/FeO ratios. An example of this is 
given by the phase diagram MgO—FeO.Fe,O0,—SiO, (fig. 9) drawn to repre- 
sent phase relationships in the system MgO—iron oxide—SiO, in air (Muan 
and Osborn, 1956). Points located exactly in the plane chosen for projection 
have Fe.0,/FeO ratios of 2.22 (on a weight “¢ basis). Mixtures ete on 
the Fe,O, side of this plane have Fe.O,/FeO ratios larger than 2.22, and 
mixtures located on the FeO side of the plane have Fe.O,/FeO ratios below 
this value. 


Similarly, approximate locations of fractionation curves may be indi- 


cated by projecting them into a chosen plane. An example of this is shown in 


figure 10, representing phase equilibria in the system FeOQ—Fe,0,—AI_0 


SiO. in air (Muan, 1957 
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Fig. 9. Diagram to illustrate phase relationships in the system MgQ-—FeO—-Fe.0, 
SiO, in air. For the sake of simplicity the relationships are pictured as a ternary system 
with MgO, FeO.FecO,, and SiO, as components, Heavy lines are boundary curves (these 
curves are not boundary curves in the strict meaning of the term as applied to ternary 
systems but rather traces of the curved boundary surfaces between two adjacent primary- 
phase volumes in the “plane” which the diagram represents), with arrows pointing in 
the directions of falling temperatures, and light solid lines are liquidus isotherms at 
100°C intervals. Boundary curves and isotherms are drawn as dashed lines in regions 
where quenching data were not obtained. Light dash-triple dot lines are curves approxi- 
mately passing through points of equal Fe.O,/FeO ratios at liquidus temperatures. Me- 
dium lines with stippling on one side indicate limit of two liquid region. 


PATHS OF EQUILIBRIUM CRYSTALLIZATION 

Phase equilibrium relationships can best be described by derivation of 
paths of equilibrium crystallization. This term has been used primarily by 
petrologists concerned with the differentiation of a magma into various min- 
erals and rocks, but because of its general usefulness in phase equilibrium 
considerations it deserves much wider attention. As originally defined the 
path of equilibrium crystallization means the course along which the liquid 
composition changes as heat is withdrawn from the system under equilibrium 
conditions. It also, however, implies the changes taking place in the composi- 
tion of the gas and crystalline phases during crystallization, In a wider mean- 
ing, therefore, the path of equilibrium crystallization is a continuous 


description of the changes in phase assemblage taking place as the tempera- 
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Fig. 10. Diagram of the system iron oxide Al,O;—SiO, in air presented in order to 
illustrate relationship of composition of liquid and crystalline phases in equilibrium. 
Heavy lines are boundary curves and medium dashed lines are fractionation curves. Lines 
with stippling on one side indicate limit of two liquid region. (Diagram after Muan 
(1957) ). 
ture is lowered from liquidus temperature, and equilibrium maintained among 
the phases. The same sequence of phase changes of course occurs upon 
heating under equilibrium conditions, only in the opposite direction, and one 
could equally well talk about a “path of melting”. Thus a knowledge of the 
path of equilibrium crystallization in a system enables one to determine the 
composition of all phases present and their relative amounts at any tempera- 
ture within the interval in which the system has been investigated. 

In most dry oxide systems, where the gas phase can be neglected, the 
path of equilibrium crystallization can be derived on the assumption that the 
total composition of condensed phases remains constant. In systems where 
changes in oxidation states take place, this relationship usually does not hold. 
In reacting with the atmosphere, as temperature is changed, oxygen is added 


to or subtracted from the condensed phases concomitant to a proportional 
change in the MeO,/MeO, ratio. The total composition of the condensed 
phases thus changes continuously during crystallization of the liquid, but only 
with respect to the relative proportions of the two components MeQ, and 
MeO,. This change can be described by straight lines passing through the O 
corner of the system in question. By combining these lines with liquidus iso- 
therms, fractionation curves, O, isobars or lines of equal pco./pr, ratios, 
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one is able to derive paths of equilibrium crystallization for various idealized 
conditions, which conditions in practice approach more or less closely. In the 
following, methods for derivation of paths of equilibrium crystallization in 
binary, ternary and quaternary systems involving changes in oxidation states 
will be reviewed briefly, using the systems Fe—O, FeOQ—Fe,0,—SiO, and 
MgO—FeO—Fe,0,—SiO, as examples. Four different idealized conditions 
will be considered: 


(1) Constant total composition of condensed phases 

(2) Constant O, pressure 

(3) Constant mixing ratio peo,/p=, of the gas phase 

(4) Crystallization in contact with metallic iron 

The first of these conditions is closely approximated in a closed system 
when the volume of the gas phase is small. Conditions under (2) and (3) are 
easily realized experimentally and have been used extensively in recent phase 
equilibrium studies (Muan, 1955; Muan and Osborn, 1956; Muan and Gee, 
1956; Muan, 1957a, b). Finally, the fourth condition obtains under the 
experimental conditions first used by Bowen and Schairer for studying iron 
silicate systems (Bowen and Schairer, 1932). Throughout the following dis- 
cussion it is assumed that perfect equilibrium exists among the phases. Also, 
it is assumed that no reaction takes place between the sample and the material 
of the crucible in which the sample is contained, except where iron crucibles 
are used (case 4 above). 

The System Fe—O.—In a binary system the situation is usually so 
simple that the term path of equilibrium crystallization is rarely used at all. 
However, for the sake of completeness and because the present system is more 
complicated than ordinary dry binary systems, the relationships will be dis- 
cussed briefly. Under conditions of constant total composition of condensed 
phases (case 1) the path of equilibrium crystallization can be read from the 
diagram in figure 4 by following a vertical (constant composition) line from 
high to low temperatures, For instance, a mixture of composition 24 weight 
“> oxygen will go through the following sequence of phase changes as heat is 
withdrawn from the closed system under equilibrium conditions. Above 
1395°C the mixture is all liquid. At 1395°C the first traces of wiistite appear, 
and wiistite and liquid coexist in equilibrium over the temperature range 1395 
to 1388°C. At the latter temperature the last trace of liquid disappears, and 
wiistite is the only condensed phase present in the interval 1388-740°C. Mag- 
netite starts forming at 740°C, and this phase coexists in equilibrium with 
wiistite down to 560°C. At this temperature, corresponding to an invariant 
situation, the four phases wiistite, magnetite, metallic iron, gas are present 
together in equilibrium. Finally, below 560°C metallic iron and magnetite 
are the stable condensed phases, During the crystallization process the O, 
pressure changes continuously in a manner which may be read directly off the 
diagram in figure 4. The intersection of the vertical (constant composition) 
line and the O, isobar at any temperature gives the equilibrium O, pressure 
of the gas phase. 


At constant O, pressure (case 2) the path of equilibrium crystallization 
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coincides with the O, isobar corresponding to the chosen O, pressure. As an 
example, consider crystallization at a constant O, pressure of 10—* atm. Above 
1626°C a metallic liquid is the only condensed phases present. At 1626° a 
metallic liquid and an oxide liquid coexist in equilibrium with the gas 
phase. Between 1626 and 1397°C liquid oxide is the only condensed phase 
existing under equilibrium conditions. At 1397°C liquid and crystals of 
wiistite and gas are present together in equilibrium. In the temperature in- 
terval 1397-1280°C wiistite is in equilibrium with the gas, and at 1280°C 
wiistite and magnetite coexist with the gas. From 1280 to 840° magnetite 
is the stable phase in an atmosphere of this O, pressure. At 840°C, magnetite 
and hematite coexist, and below 840°C hematite is in stable equilibrium with 
the gas. One important difference between this case and the case discussed 
under (1) should be emphasized: In the present case two condensed phases 
(such as for instance liquid and wiistite) can coexist in equilibrium only at 
one temperature, This is because we have “used up” one degree of freedom by 
choosing the O, pressure. Under the conditions stipulated under (1), how- 
ever, no restrictions in the sense of the phase rule are imposed on the system, 
and two condensed phases can coexist over a temperature interval. 

The third case, crystallization under conditions of constant peo,/px, 
ratios of the gas phase, is very similar to the case just discussed and needs 
no special description. Lines of equal peo,/pu, ratios are paths of equilib- 
rium crystallization under these conditions. 

If the crystallization takes place in contact with metallic iron (case 4), 
the changes in phase assemblage can be read from the diagram (fig. 4) 
by following the boundary curve between metallic iron and the lowest stable 
oxide, from high to low temperatures. Between 1535 and 1524°C solid iron 
is in equilibrium with a metallic liquid and gas. At 1524°C solid iron, a 
metallic liquid and an oxidic liquid coexist in equilbrium with gas, Between 
1524 and 1371°C liquid oxide is in equilibrium with iron and gas. In the 
temperature interval from 1371 to 560°C wiistite is the equilibrium condensed 
phase in contact with iron, and below 560°C magnetite coexists in equilibrium , 
with iron. Again, the variation in O, pressure during crystallization is deter- 
mined as the intersections between the O, isobars and the boundary curve 
along which iron is one of the equilibrium phases. 

The System FeO—Fe,0,—SiO,.—Paths of equilibrium crystallization 


under conditions of constant total composition of condensed phases (case 1) 


can be derived by methods analogous to those used for ordinary dry silicate 


systems. Such methods have been described in detail in many papers emanat- 
ing from the Geophysical Laboratory of the Carnegie Institution of Washing- 
ton.'? The application of these principles to the system FeQ—Fe,0,—Si0, 
has been described in detail elsewhere (Muan, 1955). Only one example will 
be given here for the purpose of illustrating the general relationships. The 
crystallization paths can be derived from diagrams showing liquidus isotherms 
(fig. 5), O. isobars along liquidus surface (fig. 6) and fractionation curves 


See for example Bowen's classic book (N. L. Bowen, Evolution of the igneous rocks, 
reprinted in 1956 by Dover Publications, Inc., New York, 332 p.) 
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(fig. 8) on the basis of the assumption that the total composition of con- 
densed phases remains constant. 


vases, Heavy lines are boundary curves, medium solid line equilibrium crystal- 


> 

WEIGHT PER CEN Fe,0, z3 
Fig. 1] Diagram of a part of the systemFeO—Fe,0 SiO, shown in order to illus- 
rate paths of equilibrium crystallization under conditions of constant total composition of 

ondensed pl 


izalio 


n curve and medium dash lines fractionation curves. Light solid lines are liquidus 
isotherms, light dash-dot lines O, isobars for points on liquidus surface, and light dash 
ines conjugation lines. (Diagram after Muan (1955) ). 

On the basis of these diagrams, paths of equilibrium crystallization for 
a part of the system FeOQ—Fe,0,—SiO, are shown in figure 11. Heavy lines 
are boundary curves. Medium dash lines are fractionation curves, medium 
solid lines are equilibrium crystallization curves and light dash lines are 
conjugation lines. Light solid lines are liquidus isotherms and light dash-dot 
lines are O, isobars. 

Mixture D is chosen to illustrate paths of crystallization in the field of 
wiistite. At approximately 1370°C wiistite starts separating out and the liquid 
composition moves along the equilibrium crystallization curve D-d (mediun. 
solid line). The composition of the wiistite in equilibrium with liquid, deter- 
mined by the tangent to the fractionation curve through the points represer*- 
ing liquid composition, changes continuously from d” (31 weight % Fe,O,) 
to d’” (28 weight “ Fe,O,). Fayalite starts separating out together with 
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wiistite from the liquid at point d, and the liquid composition moves along the 
boundary curve wiistite-fayalite towards d’. From the point d where the 
boundary curve is reached, a reversal occurs in the direction of change in 
composition of wiistite. the Fe,O, content increasing as the liquid composition 
moves along the boundary curve. The liquid disappears at point d’, the end 
product being a mixture of fayalite and wiistite. The composition d”” of 
wiistite at this point, determined by the conjugation line d’-d’”’, falls on the 
extension of the straight line from 2FeO.SiO, through point D. During the 
crystallization process, the O, pressure changes continuously. It first decreases 
along the curved path D-d from approximately 10~* atm. at D to 10~*°* atm. 
at d. From here on, with fayalite and wiistite coprecipitating, the change in 
O. pressure reverses itself, increasing to approximately 10~—*°-* atm. at point d’. 

In contrast to this behavior, consider the path of crystallization at con- 
stant O, pressure. This case, likewise, has been discussed in detail in a pre- 
vious publication (Muan, 1955). Under these conditions of crystallization, 
the composition of the condensed phases does not in general remain constant 
during the process. By reacting with the gas phase, oxygen is continuously 
added to or subtracted from the condensed phases as temperature varies, 
concomitant to a proportional change in Fe,O,/FeO ratio of condensed phases, 
while the amount of the other constituent (SiO,) remains unchanged. These 
compositional changes are described by straight lines pointing to the O corner 
of the triangle representing the system Fe—Si—O, of which FeOQ—Fe,0 
SiO, is a part. In the previous discussion of this system (Muan, 1955) these 
lines were referred to as iso-silica lines to designate the constant amount of 
SiO, present in the condensed phases during the reaction between condensed 
phases and gas. These lines have the same meaning as the curves referred to 
as dissociation paths by Richards and White (1954). Beacuse oxygen can 
either be added to or subtracted from the condensed phases, the term oxygen 
reaction lines seems more appropriate and will be used as a general term to 
designate these lines in this and in future papers from this laboratory. An 
example of path of crystallization at constant O, pressure is illustrated in 
figure 12. Here the same line symbols have been used as before and in addi- 
tion dash-cross lines are used to indicate changes in total composition of 
condensed phases. 

In mixture C wiistite starts separating out at about 1350°C, and as heat 
is withdrawn, the liquid composition moves along the wiistite liquidus surface 
following the 10~°* atm. isobar. The wiistite composition changes continuously 
in a way which can be derived from the fractionation curves. These curves 
are not drawn in figure 12, in order to avoid the conius’ 1 caused by ive 


many lines. When the liquid composition reaches point c, wiistite disappears 
and magnetite starts separating out with temperature remaining constant. 
When the last trace of wiistite disappears, the liquid composition starts moving 
along the magnetite liquidus surface following the 10~-* atm. isobar, with 
magnetite continuously separating out. When the boundary curve magnetite- 
tridymite is reached, at point c’, the liquid disappears, leaving as end product 
a mixture of magnetite and tridymite. The total composition during this 
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FeO Fe,0, 
WEIGHT PER CENT 
Fig. 12. Diagram of the system FeOQ—Fe.0,—SiO, shown in order to illustrate paths 
of equilibrium crystallization at constant O, pressure. In addition to the line symbols 
explained for fig. 11 light dash-cross lines are used to represent changes in total composi- 
tion of condensed phases during crystallization, and lines with stippling on one side in- 
dicate limit of two liquid region. (Diagram after Muan (1955) ). 


process changes along the straight line from C to C’. The sudden discontinu- 
ous change in total composition of condensed phases as a boundary curve is 


reached (points ¢ and c’ in fig. 12) is analgous to the change along a hori- 
zontal line in a two phase (two condensed phases) region in the system Fe— 
O (fig. 4). This situation is caused by the requirement of disappearance of 
one phase before the temperature can change, due to the restrictions imposed 
on the system resulting from our choosing the O, pressure. 

An O, isobaric section through the ternary system FeQ—Fe,0,—SiO, 
is not a true binary system, because the compositions of all phases can not be 


represented in terms of two chosen components. However, aside from this 
limitation, the system can be treated like a binary system as far as the crystal- 
lization process is concerned, The liquid compositions will always be repre- 
sented by points along the O, isobar in question, regardless of the nature of 
the crystalline phases separating out, whether they be pure phases or solid 
solution crystals. The O, of the atmosphere acts as a regulator to bring the 
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liquid composition back to its equilibrium value in spite of the crystallizing 
solids trying to bring the liquid compositions off the isobar. The lowest tem- 
perature at which a liquid can exist in equilibrium under these conditions, 
therefore, is the lowest liquidus temperature along the O, isobar correspond- 
ing to the chosen atmosphere. 

Crystallization under conditions of constant peo./p», ratio (case 3) is 
so similar to the case considered above that no special discussion of this case 
is necessary. The relationships are analogous to the above, with lines of equal 
peo, ‘ps, ratios replacing O, isobars (fig. 7). 

Crystallization in contact with metallic iron is a situation corresponding 
to the experimental conditions used by Bowen and Schairer (1932). The path 
of equilibrium crystallization in the FeOQ—Fe,0,—SiO, system under these 
conditions are the boundary curves between the field of metallic iron and, with 
increasing SiO, content, wiistite, fayalite, tridymite, cristobalite, in that order 
(figs. 5-7). A mixture whose composition is located at the boundary curve 
metallic iron—wiistite behaves as follows: Above liquidus temperatures’* a 
liquid oxide phase is in equilibrium with iron. As temperature is lowered be- 
low the liquidus, wiistite crystallizes out, with liquid composition changing 
along the boundary curve. At 1177°C the liquid is also saturated with fayalite, 
and this phase starts separating out, with temperature remaining constant as 
heat is withdrawn from the system, As the last trace of liquid disappears, 
temperature decreases, with wiistite and fayalite coexisting with metallic iron 
below 1177°C (but above 560°C). In a mixture of composition located at 
the boundary curve between the fayalite composition’* and the 1177° peri- 
tectic, fayalite crystallizes out, with liquid composition changing along the 
fayalite-iron boundary curve until the temperature of 1177°C is reached. At 
this temperature the liquid disappears, leaving as end product a mixture of 
fayalite and wiistite in equilibrium with metallic iron. A mixture with a com- 
position corresponding to fayalite'* is all liquid above 1205°C (Bowen and 
Schairer, 1932). At this temperature the mixture crystallizes completely to 


fayalite in equilibrium with metallic iron. Similarly, any mixture of composi- 
, 


tion on the SiO, side of the fayalite composition crystallizes to a mixture of 
fayalite and tridymite at the peritectic temperature of 1178°C. Fayalite is the 
primary crystalline phase in mixtures located between 2FeO.SiO, and this 
peritectic, while tridymite or cristobalite is the primary phase in mixtures 
located between SiO. and the peritectic composition. The Ss) stem iron oxide 

SiO, in contact with metallic iron is not truly binary, because compositions 
of condensed phases can not be represented in terms of two chosen com- 
ponents. Aside from this limitation, however, the path of crystallization can 
be treated in a manner analogous to that for a binary system. The usual 
representation of the system in terms of a projection into the join FeOQ—Si0, 
as illustrated in figure 13, after Bowen and Schairer (1932), is therefore very 
instructive and quite sufficient for most purposes. It should be pointed out 


Liquidus temperature is used here to indicate highest temperature of stable existence 
f a crystalline oxide phase 


The exact location of this composition is determined as the point of intersection be- 
tween the boundary curve and the oxygen reaction line passing through the 2FeO.SiO, 
composition point. 
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specifically, however, that the apparent eutectics in the join FeOQ—SiO, are in 
reality peritectics in the ternary system FeQ—Fe,0,—Si0,. 


CRISTOBAL ITE + LIQUID 


TEMPERATURE 


FAYALITE + TRIDYMITE 


| 
FeO 2FeO SiO, 


x 80 
wt 

Fig. 13. Diagram illustrating phase equilibrium relationships at liquidus tempera- 
tures in the system iron oxide—SiO, in contact with metallic iron. The diagram has been 
plotted with the appearance of a binary system with FeO and SiO, as components. (After 
Bowen and Schairer (1932)). 

The System MgO—F eO—F e,G,—SiO,.— Crystallization under conditions 
of constant total composition of condensed phases is analogous to the situa- 
tion in ordinary dry oxide systems. Methods for describing path of crystalliza- 
tion in dry quaternary silicate systems have been discussed in the literature.*® 
Applications of such methods to the system MgOQ—FeO—Fe,0,—SiO, were 
discussed in detail in a previous paper from this laboratory (Muan and 
Osborn, 1956). Suffice it here to summarize the relationships briefly: In 
figure 14 the system MgO—FeO-—Fe,0,—SiO, is illustrated as a regular 
tetrahedron. Quaternary univariant lines along which three crystalline phases, 
one liquid and a gas phase coexist in equilibrium, extend through the tetra- 
hedron from the left hand face to the right hand face. As an example of 
crystallization at constant total composition consider a mixture of composi- 


Compare footnote 12. 
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3S 


Fig. 14. Diagram illustrating phase relationships in the system MgO—FeO—Fe.0, 
SiOz, after Muan and Osborn (1956). The system is pictured as a regular tetrahedron 
with one component at each apex. Compounds, with compositions given in terms of the 
components, are plotted in weight per cent. The plane MgO—FeO—Fe,0,—Si0, is indi- 
cated with light solid lines. Boundary curves are drawn as medium solid lines in the base 
(the system MgO—FeO—SiO.) and in the right-hand face (the system FeO—Fe,O,; 
SiO.) and as light lines in the plane Mg0—FeO.Fe,0,—SiO, (corresponding approxi- 
mately to the system MgOQ—FeO—Fe,0,—SiO, in air). Quaternary univariant lines, ex- 
tending into the tetrahedron from the bounding ternary systems, are also drawn as medium 
solid lines. Medium dashed lines represent inferred boundary curves in the bounding 
ternary systems MgO—Fe.0,—SiO, and Mg0—FeO—Fe.Os. Medium lines with stippling 
on one side indicate limit of region of coexistence of two silicate liquids. 


tion ¢ as illustrated in figure 15, This point is located on the univariant line 
along which magnesioferrite, pyroxene, olivine, liquid and gas exist together 
in equilibrium. Immediately above liquidus temperature a liquid of composi- 
tion c is in equilibrium with gas. At liquidus temperature, magnesioferrite 
(composition c’), pyroxene (composition c”) and olivine (composition c””), 
start separating out, and as heat is continuously withdrawn from the system, 
the liquid composition changes to the right along the boundary curve from c 
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toward m. As crystallization proceeds, the three crystalline phases change 
compositions along the joins MgO.Fe,0,-—FeO.Fe,0,, MgO.Si0,—Fe0.Si0,, 
and 2Mg0.Si0,—2FeO0.SiO, respectively, and the O, pressure decreases con- 
tinuously. When the liquid composition reaches point m, the crystalline phases 
in equilibrium have the compositions m’, m”, and m’” respectively. These 
points are located so that point c, representing the constant total composition, 
falls in the plane m’-m”-m’”’. Complete crystallization therefore takes place 
when the liquid reaches point m. The final product is a mixture of magnesio- 
ferrite of composition m’, pyroxene of composition m”, and olivine of com- 
position m’””. 


Fig. 15. Diagram to illustrate paths of equilibrium crystallization for mixtures in 
the system MgO—FeO—Fe,0,—SiO, under conditions of constant total composition of 
condensed phases. Open circles and solid dots represent compositions of liquid and 
crystalline phases, respectively, for mixtures whose paths of crystallization are discussed 
in the text. Light solid lines indicate the joins MgO.Fe,0;—FeO.Fe,0,, MgO0.Si0,— 
FeO.SiO., and 2MgO.SiO,—2FeO.SiOs. The two shaded triangular planes are formed by 
lines connecting points representing compositions of crystalline phases in equilibrium 
with liquid at the quaternary univariant situation. (Diagram after Muan and Osborn 
(1956)). 


A mixture located in the plane g’-g”-g”” will crystallize completely only 
after the liquid composition has reached point g. Here g’, g” and g”” represent 
compositions of magnesioferrite, pyroxene and olivine, respectively, that to- 
gether with tridymite are in equilibrium with melt of composition g at the 
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quaternary liquidus invariant point. Thus, under conditions of constant total 
composition of condensed phases, the lowest temperature of liquid existence 
in a quaternary system is the temperature of a liquidus invariant point 
(eutectic or peritectic) or the temperature of that point on a quaternary 
univariant line representing composition of the liquid phase in equilibrium 
with crystalline phases whose composition points form a plane in which the 
composition of the original starting mixture is located. 


Fe,0, 


a 
FeO Fe,0, 


MgO 
Fig. 16. Diagram to illustrate paths of equilibrium crystallization for mixtures in the 
system MgO—FeO—Fe,0;—SiO, in air, Open circles and solid dots represent composi- 


tions of liquid and crystalline phases, respectively. Light solid lines indicate the joins 
Mg0.Fe.0;—FeO.Fe.0s, and 2MgO.SiO,—2FeO.SiO, and outlines 
of the plane MgO—FeO.Fe.0;—SiO:. Medium solid lines are boundary curves. The shaded 
planes are formed by lines connecting points representing compositions of crystalline 
phases. (Diagram after Muan and Osborn (1956) ). 


In contrast to this behavior, consider crystallization at constant O, pres- 
sure, as demonstrated by a specific case illustrated in figure 16. Crystallization 
in air is chosen as an example, i.e. at a constant O, pressure of 0.21 atm. 
Compositions of mixtures at liquidus temperatures in air are located in an 
irregularly curved surface, the 0.21 atm. O, isobaric surface, which extends 
into the tetrahedron from the corresponding isobaric line in the system FeO 
Fe,0,—SiO, (fig. 6). For the sake of simplicity this O, isobaric surface has 
been pictured as a plane MgO—FeO.Fe,0,—SiO, in figure 16. Consider a 
mixture of composition represented by point N on the boundary curve between 
tridymite and magnesioferrite in this plane. At the liquidus temperature for 
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this mixture, liquid of composition N is in equilibrium with infinitesimal 
amounts of tridymite (SiO,) and magnesioferrite of composition n’. With 
magnesioferrite of changing composition and tridymite crystallizing out to- 
gether, the liquid composition moves along the boundary curve in the surface 
MgO—Fe0.Fe,0,—SiO,. The total composition during this process changes 
along the oxygen reaction line’* through the point of original composition 
from N in the direction toward N’. Complete crystallization takes place when 
the liquid composition reaches point p, at which point the magnesioferrite 
has a composition p’. This latter point is located so that a line from the silica 
corner to p’ goes through N’, the point of intersection of the oxygen reaction 
line with the plane MgO.Fe,0,—FeO.Fe,0,—SiO,. The final product is thus 
a mixture of tridymite and magnesioferrite of composition p’, the relative 
amounts being determined by the ratio N’p’, SiO,N’. 

Similarly it can be shown that all mixtures in the plane Mg0—FeO.Fe,0, 

SiO, located on the SiO,—FeO.Fe,0, side of the dashed curve from SiO, 
to q will crystallize completely to a mixture of magnesioferrite and tridymite. 
The significance of the dashed line SiO,-q is the following: As the liquid com- 
position changes along the tridymite-magnesioferrite boundary curve from 
r to d, the magnesioferrite composition changes continuously from FeO.Fe,0, 
to d’ (d’ is the composition of magnesioferrite in equilibrium with liquid d). 
For mixtures in the Mg0—FeO.Fe,0,—SiO, plane to crystallize completely 
to magnesioferrite and tridymite, the oxygen reaction line through the point 
of starting composition must intersect the plane SiO,—FeO.Fe,0,—d’. In 
other words, the oxygen reaction lines through the SiO, d’ edge of this plane, 
when cutting through the plane Mg0—FeO0.Fe,0,—SiO, along the line SiO,- 
q, separate the mixtures into two groups with respect to crystallization paths. 
Mixtures located on the SiO,—FeO.Fe.0, side of the line SiO,-q crystallize 
completely to a mixture of magnesioferrite and tridymite at a temperature 
corresponding to the liquidus for a mixture of composition located on the 
boundary curve r—d. Mixtures located on the other side of the line SiO,-q 
have oxygen reaction lines that intersect the plane SiO,—d’—d”, d” being 
the composition of the pyroxene in equilibrium with liquid of composition d. 
The final product for these mixtures, therefore, is magnesioferrite, tridymite 
and pyroxene, and the last trace of liquid disappears at the liquidus tempera- 
ture for a mixture of composition d. 

Noteworthy among the relationships explained above is the fact that the 
liquid composition during crystallization remains in the 0.21 atm. O, surface 
(here pictured as the plane MgO—FeO.Fe,0,—SiO,). The reason for this is 
the same as explained for FeOQ—Fe,0,—SiO,. The O, of the atmosphere com- 
pensates for the tendency of the crystallizing solid phases to push the liquid 
composition out of the isobaric surface. Hence, at constant O, pressure, a 
liquid phase can never exist in equilibrium at a temperature lower than the 
lowest liquidus temperature in the O, isobaric surface corresponding to the 


chosen atmospheric conditions, This temperature is usually that correspond- 


‘* The oxygen reaction lines have a direction pointing towards the O apex of the tetra- 
hedron representing the system Mg—Fe—Si—O, and have the same significance as the 
corresponding lines for the system FeOQ—Fe,0;—SiO, (compare page 194). 
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ing to the intersection between the O. isobaric surface and a quaternary 
univariant line. Only on the O, isobaric surface which passes through the point 
representing composition of the liquid at a quaternary liquidus invariant 
situation will a liquid phase exist in equilibrium with four crystalline phases 
at the temperature of this situation. The isobaric surfaces are not true ternary 
systems, because the compositions of the phases can not be represented by 
three components. However, in spite of this limitation, in the derivation of 
paths of equilibrium crystallization at constant O, pressure, they can be 
treated like ternary systems. 

The crystallization process at constant CO, /H, ratios requires no detailed 
discussion. The relationships are analogous to above, with lines of equal 
CO./H, ratios replacing the O, isobars. 

The paths of equilibrium crystallization of mixtures in the system MgO 
FeO—Fe,0,—SiO, in contact with metallic iron can be derived by methods 
analogous to those described for the systems Fe—O and FeQ—Fe,0,—Si0O.. 
The small primary field of metallic iron in the FeQ—Fe,0,—SiO, system 
(figs. 5-7) extends as a thin phase volume into the tetrahedron MgO—FeO— 
Fe.0,—SiO,, adjacent to the basal plane MgQ0—FeO—SiO,. This phase 
volume was left out of figures 14-16 in order to simplify the drawings. A pure 
iron oxide liquid in contact with metallic iron contains approximately 10 
weight % Fe.O,"" at liquidus temperature. As SiO, is added, the FeO, /FeO 
ratio of the liquid in equilibrium with iron (see boundary curve between 
metallic iron and oxide phases in figs. 5-7) decreases. The decrease is caused 
by the increase in acidity'* of the liquid as SiO, content increases. This effect 
more than balances the opposite effect on Fe,O,/FeO ratio of the decrease in 
liquidus temperatures caused by the first SiO, additions. As more SiO, is 
added, liquidus temperatures vary relatively little with SiO, content in the 
region where fayalite (2FeO.SiO.) is the primary crystalline phase. When the 
liquid becomes saturated with SiO,, liquidus temperatures rise sharply as 
SiO, content increases, This increase in liquidus temperatures magnifies the 
acidity effect of SiO, on FeO, FeO ratios, and hence the boundary curve 
between metallic iron and cristobalite approaches the join FeQ—SiO, as the 
temperature of melting of iron is reached. As MgO is added to silica-iron 
oxide mixtures, the basic character of MgO may be expected to increase the 
Fe.0,/FeO ratio of the liquid phase. However, the MgO addition at the same 
time raises liquidus temperatures very markedly. The Fe,0, content of 
liquids in equilibrium with metallic iron is therefore very low for all composi- 
tions in the system MgO—iron oxide—SiO, except those close to the iron 
oxide apex of the tetrahedron representing the system. The relationships can 
be seen quantitatively from the analytical data of Bowen and Schairer (1935). 
The FeO, content is less than 2 weight “ for most compositions, and hence 
** According to Darken and Gurry (1946) the FeO, content is 8.22 weight %, while 
Bowen and Schairer (1932) give a value of 11.56 weight %. 


* A quantitative acid-base scale is not available for melts like we have for aqueous solu- 
tions. The term acid is used here to describe the behavior of an ion of high field strength 
(Si), in accordance with modern concepts of acid-base relationships in oxide systems 
(Weyl, 1956). 
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Fig. 17. Diagram illustrating phase relationships in the system MgO—FeO—Si0O,, 
after Bowen and Schairer (1935). Heavy lines are boundary curves with arrows pointing 
in the directions of falling temperatures. Light solid lines are liquidus isotherms at 
100°C. intervals. Boundary curves and liquidus isotherms are dashed in regions where 
liquidus data were not obtained. Medium lines with stippling on one side indicate limit 
of two-liquid region. 
the iron phase volume extends only as a very thin slice into the tetrahedron 
of figures 14-16, adjacent to the basal plane. 

It will be recalled that in the system FeOQ—Fe,0,—SiO, the path of 
crystallization for mixtures in contact with metallic iron was the boundary 
curve between metallic iron and the adjacent phase areas. In the present case, 
similarly, liquid compositions during crystallization must remain in the 
boundary surface between the primary phase volume of metallic iron and the 
adjacent phase volumes in the quaternary system, magnesio-wiistite, olivine, 
pyroxene, tridymite, cristobalite, depending on composition of the mixture. 
Projections of these phase volumes are indicated in figure 17, which is a 
reproduction of the system MgO—“FeO”—SiO, after Bowen and Schairer 
(1935). Their diagram is obtained by projecting actual compositions of 


mixtures located along the boundary surface limiting the iron phase volume, 
into the base (MgQ—FeO—SiO,) of the tetrahedron representing the MgO— 
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FeO—Fe,0,—SiO, system. Analogous to the situation in an O, isobaric sec- 
tion through the quaternary system, liquid compositions in the present case 
remain in a fixed surface, namely the surface bounding the phase volume of 
metallic iron. While in the former case reaction with oxygen of the gas phase 
was the mechanism for keeping the liquid composition in the appropriate 
surface, the same is accomplished in the present case by reaction of the oxide 
melt with the iron of the crucible in which the mixture is contained. Regard- 
less of the nature of the crystalline phases separating out, the liquid composi- 
tion remains in the surface bounding the phase volume of metallic iron, The 
system is not truly ternary under these conditions, because liquid compositions 
in the irregularly curved surface can not be represented in terms of three 
chosen components. However, aside from this limitation it can be treated like 
a ternary system. Because of the closeness of the iron boundary surface to the 
basal plane, Mg0—FeO—Si0O,, the paths of crystallization can be discussed 
very satisfactorily in terms of a projection into this plane. The diagram ob- 
tained in this way has the appearance of a ternary system, and was treated 
as such by Bowen and Schairer. The reader is referred to their paper for a 
detailed discussion of paths of equilibrium crystallization. It should be kept 
clearly in mind, however, that the condensed phases contain some Fe**, and 
that metallic iron is present as a phase in addition to the oxide phases, Hence 
apparent ternary invariant situations in the projected diagram are actually 
quaternary liquidus invariant situations with metallic iron as one of the 
crystalline phases present. 


SUMMARY 


Phase equilibrium studies of oxide systems at high temperatures have 
been restricted in the past largely to systems containing only oxides of noble 
gas type ions such as for instance CaO, MgO, Al,O,, SiO,. Data are essential- 
ly lacking on oxide systems involving elements where changes in oxidation 
states occur. The study of the latter systems is complicated because the com- 
position of the gas phase must be controlled carefully and varied between 
wide limits. At reasonably high levels of O, pressures (1 atm.—10~* atm.) 
such equilibria may be studied experimentally in O, atmospheres where the 
desired pressure is adjusted by vacuum techniques or by mixing O, with an 
inert gas. When the partial pressure of O, becomes low (<10~-* atm.) the 
gas phase is so diluted in O, that a very substantial volume of gas would be 
necessary to change the composition of the condensed phases to the equilib- 
rium value. With a moderate rate of flow of gas through the furnace in which 
the equilibrium study is made, it is therefore not possible to reach equilib- 
rium within a reasonable period of time. In such cases indirect methods of 
obtaining the desired O, pressure must be used. The most convenient method 
is to use a mixed gas system containing an O, bearing species which upon 
heating dissociates under liberation of O,. Convenient gases for this purpose 
are CO, or H,O. The degree of dissociation and hence the O, pressure may 
be controlled by mixing CO or H, with the gas in controlled proportions. 

Most of the few available phase diagrams for oxide systems involving 
changes in oxidation states are very incomplete. In addition to the usual 
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representation of composition of condensed phases such diagrams should 
contain information regarding the equilibrium composition of the gas phase 
as well. It is suggested that binary systems be represented in an usual com- 
position-temperature diagram with O, isobars or lines of equal CO,/H, equilib- 
rium ratios superimposed to indicate gas composition. In ternary systems 
showing phase relationships at liquidus temperatures the ordinary projection 
of liquidus surface with liquidus isotherms should be supplemented by di- 
agrams showing O, isobars (or lines of equal pco./px, ratios) along the 
liquidus surface and fractionation curves. Quaternary systems are difficult to 
represent graphically, Phase relationships at liquidus temperatures are best 
illustrated by projecting compositions from their actual locations in irregular- 
ly curved O, isobaric surfaces into chosen planes. In the resulting diagrams, 
which have the appearance of ternary systems, the actual compositions of 
mixtures are indicated by lines of equal ratios (MeO,/MeO,) of the oxides 
of the element with changing oxidation state (such as for instance Fe,0,/ 
FeO). 

Path of equilibrium crystallization is a term used for the continuous 
description of changes in phase assemblage taking place as heat is withdrawn 
from a system and equilibrium is maintained among the phases during the 
process. Iron oxides and iron oxide containing silicates are chosen as examples 
of paths of equilibrium crystallization in systems where changes in oxidation 
states occur. Four idealized conditions of crystallization are considered, each 
one corresponding to a situation which can be closely approximated experi- 
mentally. The first condition is that of crystallization at constant total com- 
position of condensed phases. This situation is approximated in a closed 
system when the volume of the gas phase is small. The path of crystallization 
can be derived by well known methods by the use of diagrams showing liqui- 
dus isotherms, fractionation curves and O, isobars, The second condition is 
crystallization at constant O, pressure. In this case the criteria used for 
derivation of crystallization path in ordinary dry silicate systems are not ap- 
plicable, because total composition of condensed phases does not in general 
remain constant. By reaction with the gas phase oxygen is added to or sub- 
tracted from the condensed phases concomitant to a proportional change in 
the oxide ratio MeO,/MeO,. This composition change is described by oxygen 
reaction lines, which are straight lines pointing towards the O apex of the 
model representing the system in question. By choosing the O, pressure, we 
have “used up” one degree of freedom, and the system possesses one degree 
of freedom less than in the first case considered, for the same number of 
phases present. In a binary system, therefore, two condensed phases, such as 
for instance one crystalline and one liquid phase, can coexist only at one 
temperature at a chosen O, pressure. In a ternary system the O, isobars along 
the liquidus surface are the paths along which the liquid composition changes 
during crystallization, The lowest temperature of stable existence of a liquid 
phase is usually the intersection between a boundary curve and the O, isobar 
in question. In a quaternary system, similarly, the liquid composition remains 
in the particular O, isobaric liquidus surface corresponding to the chosen O, 
pressure. The lowest temperature of liquid formation under equilibrium con- 
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ditions is usually the intersection between a quaternary liquidus univariant 
line and the O, isobaric surface in question. 

Crystallization under conditions of constant peo,/ps, ratio is analogous 
to the case above, with peo./px, ratios replacing the O, pressures in the 
discussion. 

When crystallization takes place in contact with metallic iron, the liquid 
composition is forced to move along the boundary surface (or boundary 
curve) between the phase volume (area) of metallic iron and the adjacent 
oxide phase volumes (areas). In a binary system liquid oxide and a crystal- 
line oxide phase can coexist in equilibrium with iron at one temperature only. 
In a ternary system the lowest temperature of stable existence of an oxide 
liquid is usually the temperature of a liquidus invariant situation where metal- 
lic iron is one of the phases present. In a quaternary system the lowest tem- 
perature of liquid formation, likewise, is usually the temperature of one of 
the quaternary liquidus invariant situations where metallic iron is one of the 
phases present. 
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GROUP OPERATOR VARIANCE IN 
MORPHOMETRIC WORK WITH MAPS 


RICHARD J. CHORLEY 


ABSTRACT. Two groups of inexperienced operators delimited and measured certain 
linear, areal and slope aspects of the morphometry of seven drainage basins, exhibited 
on topographic maps on a scale of 1:25,000. A statistical comparison of the combined 
group results shows that there is no significant difference between the mean values of the 
lengths, areas and channel gradients of second and third order streams obtained by each 
group, nor between the mean drainage densities. The only important source of individual 
and group variation which appeared lay in the recognition of first order stream channels. 


The application of statistical methods to the study of the geometry of 
landscape (Strahler, 1954) has made it necessary that relatively large samples 
of the values of particular morphometric units (e.g. stream lengths, basin 
areas, maximum valley-side slopes, etc.) be amssed before regions can be 
quantitatively described and before certain inter and intra-regional relation- 
ships can be established. The number of variates in each sample must often 
run into the hundreds, and this has created a purely mechanical sampling 
problem both for field workers and for those who carry out morphometric 
work with maps or aerial photographs. In the latter connection, it is interest- 
ing to determine the reliability of measurements made by operators with a 
minimum of experience and whether different operators obtain similar results 
in measuring identical landscape features. This paper concerns a comparison 
of the work of two groups of inexperienced operators, the members of which 
measured similar morphometric features from the same maps of identical 
drainage basins. 

Fourteen students, including seniors and graduate students, were given 
class instruction in the delimitation and measurement of landscape features 
from topographic maps, together with less than one hour of training in the 
operation of the chartometer and compensating polar planimeter. They were 
then allowed to divide themselves into two groups of seven students each, 
meeting at different times, and each group was assigned the maps of seven 
different drainage basins with which to work. It was purely by chance that 
the two graduate students and, generally, a majority of the better students in 
geomorphology elected to attend at the same time, forming Group 2. Each 
operator within each group was therefore concerned with a different drainage 
basin, but the two groups were each working with the same seven basins. 
Each operator worked independently and, more important, the two groups 
were physically separated. 

The drainage basins used were six fourth order and one third order, all 
located on the Dartmoor granite of Devon County, England, and depicted on 
Ordnance Survey 1:25,000 maps (Nos. 20/ 57, 58, 66, 67, 68, 76, 77 and 78). 

Apart from a minimum of initial individual assistance in drawing drain- 
age lines and watershed divides, each operator was left to delimit the 
morphometric features as he saw fit and to develop his own standards of 
accuracy of measurement. The principles of recognising first, second, third 
and fourth order stream segments and basin areas were based on those 
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modifications of Horton’s method (1945, p. 281-283) established by Strahler 
(1952, p. 1120), as illustrated by Miller (1953, fig. 12), In this revised 
method, fingertip tributaries represent first order streams and a junction of 


two streams of the same order begins a stream segment of the next higher 
order. When the measurements were complete, the final results of each group 
were combined, and comparisons of the data collected by each group were 
made. 


UNLIKE 


Fig. 1. The most unlike and like interpretations of stream channel networks by mem- 
bers of the two groups. 


Table 1 shows the total number of stream segments of each order recog- 
nised by each group. The greatest discrepancy between the discrimination of 
the two groups is to be found in their recognition of finger-tip (first order) 
tributaries. The delimitation of first order streams involves a most subjective 
element of determining what constitutes a significant crenulation in a contour 
line, thereby determining the location and extent of a headwater gully. All 
operators continued their finger-tip tributaries headward to the highest con- 
tour which, in their individual estimation, contained a significant, aligned 
crenulation. Despite variation in the interpretation of first order stream seg- 
ments, with Group 2 being apparently more assiduous in the search for head- 
water gullies (fig. 1), it is instructive to note that the total number of second 
order stream segments recognized by each group was identical, and that there 
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was only a difference of one in third orders. This similarity of results between 
the two groups can be tested statistically by applying an analysis of covariance 
(Walker and Lev, 1953) to determine whether the two sets of numbers of 
stream segments are best fitted by two separate regressions or by a single 
common line (fig. 2). In this instance the F ratio (F= .37<F 95 7.71) 
indicates that the hypothesis of a common regression line, with a regression 
coefficient of — .633, should be accepted. 


Tasre | 
Total Number of Stream Segments of each Order 


N Nz 


Group l. 405 107 
Group 2. 508 107 


Mean lengths of stream segments of each order can be calculated in two 
different ways. Firstly, each operator can accumulate on the chartometer the 
total length in inches of streams of a given order, then this figure can be 
translated into miles and divided by the number of stream segments of that 
order which were recognised. The second method involves the more arduous 
task of measuring each segment, individually changing it into miles, sum- 
ming the stream lengths and averaging them, Although the second method is 
to be preferred statistically because it provides a large sample of individual 
variates which can be treated by the standard methods of group statistics, - 
there are six sources from whence error or subjective approximation may 
enter: 

a. Setting the chartometer to zero. 

b. Measuring the length of a stream segment. 

c. Reading the measurement from the instrument. 

d. Translating the length in inches into miles. 

e. Summing the lengths. 

f. Averaging the lengths. 
In the first method: except for the initial zero setting, the single final instru- 
ment reading and the single calculation of total length inches into miles; 
sources a, c, d, and e are virtually eliminated, and for this reason it would 


seem that this method should produce the more accurate and consistent mean 
stream lengths. 


Table 2 indicates the total stream lengths obtained by the first method, 
and these have been averaged to give the mean stream lengths which are 
listed in table 3. 

Tasie 2 
Total Length of Stream Segments of each Order, in Miles 


Group 1. 111.90 45.08 35.97 
Group 2. 129.93 47.31 35.39 
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1000 


© group 1 
* group 2 


Number of Streams (N,,) 


2 3 a 
Stream Order (w) 
Fig. 2. Numbers of stream segments recognized by each group used to illustrate the 
Law of Stream Numbers, 
Taste 3 
Mean Length of Different Order Stream Segments, in Miles 


First Method. Second Method 


Order. Group 1. Group 2. Group 1. Group 2. 


.276 256 — 

421 442 400 
1.713 1.609 1.749 
3.133 3.048 3.123 


From the first, and in any case possibly the more accurate, method it is 
apparent that the mean length of first order stream segments shown on this 
scale of map of the Dartmoor region is about a quarter of a mile. The con- 
ventional chartometers which are economically available for class use have, 
as their finest graduation, divisions of one-half inch (representing, on this 
scale, approximately .200 miles) ; and are difficult to read with consistent ac- 
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curacy to within closer than .1 inch (representing about one sixth of the mean 
length of first order streams), particularly because of the small amount of 
‘play’ which is unavoidably present in the measuring wheel. It seems, there- 
fore, that the measuring of individual lengths of first order streams, the mean 
length of which is small enough to be only comparable with the smallest 
graduation on the measuring instrument, would introduce into the morpho- 
metric work just that subjective element which it is desirable to eliminate. 
The measurement of individual lengths and the calculation of means by the 
second method was limited therefore to the second, third and fourth orders, 
most of which are long enough so that instrumental and dial-reading errors 
lve only a small percentage of their length, and which, in the first 
two instances, are numerous enough to admit confidently the application of 
group statistical methods 


could inv oO 


GROUP 2 


puzZ 
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The results of these individual measurements by each group were com- 
bined (fig. 3) and strikingly exhibited the right skewness noted by Miller 
(1953) as characteristic of the frequency distributions of stream lengths and 
basin areas. As he suggested, this skewness was corrected by plotting the 
logarithms of stream lengths (fig. 4)—actually of one hundred times each 
individual length to avoid the use of negative logarithms. In order to test the 
normality of these logarithmic distributions, a chi-square test (Croxton and 
Cowden, 1939) was run on the distribution which seemed least obviously 
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normally distributed (the leagths of third order streams measured by Group 
1), employing the conventional level of significance (P) of .05 (5%). This 
test resulted in a value of y* of 2.19, indicating that, within the three degrees 
of freedom used, a normal curve is a good description of the frequency dis- 
tribution (.50<P<.70). This precaution was necessary because, in order to 
employ simple tests of similarity between mean values, it must be shown that 
the samples are representative of normally distributed populations and that 
the variances—here represented by the standard deviations (S)—of the 
samples, the means of which are to be compared, are similar. The above chi- 
square test and visual inspection of the corresponding standard deviations of 
GROUP 1 GROUP 2 
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Fig. 4. Histograms of the logarithmic values of second and third order stream 
lengths for each group. 


PYE 


the two groups (fig. 4) indicate that both these conditions are satisfied and, 
again using a level of significance of .05, ‘t’ tests were made comparing the 
means of the logarithmic values of second and third order stream lengths be- 
tween Groups 1 and 2. The statistic ‘t’, the calculation of which involves the 
arithmetic means, numbers of variates and a pooled standard deviation of 
two samples, enables one to determine the probability (P) that the observed 
difference between the two arithmetic means in question, either being the 
larger, with repeated sampling is due to chance sampling alone and does not 
represent a significant difference between the two populations from which the 
samples were drawn. Comparing, in this way, the means of the logarithmic 
values of second order (P>.90) and third order (.50<P<.60) stream 
lengths between the two groups, it was apparent—employing a critical level 
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of significance of .05—that there is no significant difference between the 
means of the logarithmic values of stream lengths of corresponding order 
between Groups 1 and 2 (fig. 4). 
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Fig. 5. Histograms of second and third order basin areas for each group (in sq. 
miles). 


Areas of first order basins in landscapes of this texture of dissection are 
difficult to obtain individually on maps of this scale. This is because of their 
difficulty of delimitation on the map, either because they are so small or be- 
cause the drawing of watershed lines on almost plane slopes becomes a highly 
subjective matter. Each group, however, delimited and measured individually 
the areas of second, third and fourth order basins; the composite mean values 
appearing in table 4 and frequency distributions for the second and third 
orders in figure 5. These distributions are strongly right-skewed and, as with 
the stream lengths, normal distribution were obtained by using the logarithms 
of the areas (fig. 6). A chi-square test on the logarithms of third order areas 
obtained by Group 2 showed them to be normally distributed (.20<P<.30). 
‘t’ tests revealed no significant difference between the means of the logarithmic 
distributions obtained by each group for either the second order basins 
(P>.90) or those of third order (P>.90). 
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Fig. 6. Histograms of the logarithmic values of second and third order basin areas 
for each group. 
Tasie 4 
Mean Values of Basin Areas, in Square Miles 


Group 1. ‘ 
Group 2. 371 


The most important morphometric measure involving both linear and 
areal aspects is the drainage density. For the region under investigation this 
was calculated in two ways (table 5). The first method was to take the mean 
of the seven calculated drainage densities for each group of seven individual 
drainage basins; the second, to divide the total measured length of channels 
in all seven basins by the sum of the seven basin areas, for each group. 
Naturally, whichever method was employed, the greater number of first order 
streams recognised by Group 2 resulted in a larger drainage density calculated 
by that group. A ‘t test involving small samples (Croxton and Cowden, 1939, 
p. 329), however, showed no significant difference (.30<P<.40) between 
the means of the drainage densities calculated by the first method. 

Slope aspects of topography are generally the most difficult and refined 
of morphometric measures to obtain from maps, involving as they do both 
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Tasie 5 
Drainage Density 


First Method. Second Method. 


Group 1. 3.24 3.17 


Group 2. 3.52 3.45 


linear and altitude considerations; and they are most subject to error, in- 
accuracy and approximation on the part of the operator. Often one must 
estimate an elevation between contour lines and sometimes linear distances 
are so short that the possible percentage error of linear measurement becomes 
quite large. The two groups obtained the gradient of each second, third and 
fourth order stream by dividing the drop in feet between each end by the 
length of the channel in miles. The mean slopes, in feet per mile, of streams 
of second, third and fourth order are given in table 6. The frequency distribu- 
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Fig. 7. Histograms of second and third order channel mean slopes for each group 
(in feet per mile). 


tions of these slopes (fig. 7) show a pronounced right-skewness, which con- 
trasts with the normal distribution of valley-side slopes (Strahler, 1950). 
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This was corrected in the usual way (fig. 8) and ‘t’ tests between the logarith- 
mic distributions obtained by the two groups showed no significant difference 
at the 5 percent level between the mean slopes of either the second order 
stream channels (P>.90) or those of third order (.20<P<.30). 


TABLE 6 
Mean Slopes, in Feet per Mile, of Stream Segments of Different Order 


Group 1. 315 
Group 2. 332 


The conclusion which one would seem to be able to draw from the above 
preliminary investigation of operator variance in morphometric work is that, 
using maps on a scale of 1:25,000 of drainage basins having a texture of dis- 
section represented by a drainage density of about 3.50, two groups of 
relatively unskilled operators collectively showed no significant difference 
between their group delimitation and measurement of the mean values of 
linear, areal and channel gradient aspects of landscape geometry—except 
where the recognition of first order stream segments was concerned. 
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mean slopes for each group. 
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DISCUSSION 
THE NOMENCLATURE OF EVAPORITE TEXTURES 
F. H. STEWART 


In a recent paper, Dr. J. T. Greensmith (1957) has snggested a number 
of new names for evaporite grains and textures, for which I feel that there is 
no real justification. 

The border between diagenesis and metamorphism is a diffuse one and is 
especially so in salts which are highly susceptible to only moderate changes 
in temperature and pressure. I have not found, however, that this raises any 
real practical difficulties in nomenclature. I do not see any objection to using 
some metamorphic textural terms in the description of evaporites (whether 
stratified or unstratified is immaterial) which have suffered widespread re- 
crystallisation or mineralogical reconstitution in response to increased tem- 
perature or pressure due to load. Evaporites commonly show excellent 
examples of changes controlled by rise of temperature on burial, classed un- 
der the well-established heading of “Geothermal Metamorphism”. The fact 
that such deposits may lie between relatively unaltered beds of limestone is 
irrelevant to their metamorphic status. 

I do not see any need for the special terms ‘evapocrysts’ and ‘neocrysts’ 
to describe individual primary or secondary evaporite grains; surely primary 
and secondary grains would do. A statement as to whether primary or second- 
ary linear or platy parallelism of grains is present or absent is better than 
using the terms “evapocrystic’, ‘evapolensic’, ‘neocrystic’ and ‘neolensic’. The 
word ‘porphyroblastic’ can be applied to metamorphosed evaporites. Would 
it really cause confusion if it were used for those in which large secondary 
crystals have formed during diagenesis? We should not try to establish a 
sharp division between sedimentary and metamorphic rocks where none 
exists, and we must avoid adding more than absolutely necessary to our 
geological jargon. 
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It is gratifying that the current method of describing evaporites should 
be so vigorously defended, as only by these means can the complexities of the 
subject be illuminated, Most petrologists would heartily agree about the 
necessity of avoiding new, and possibly cumbersome, terms. On the other 
hand, from the teaching (and students’) point of view it can be rather con- 
fusing if terms having a metamorphic connotation, such as ‘porphyroblastic’, 
and others an igneous connotation, such as ‘poikilitic’ (Dunham, 1948, p. 
220) and ‘pseudoporphyritic’ (Pettijohn, 1949, p. 355), are used to describe 
similar rocks. It is partly this inconstancy in usage which first suggested the 
need for a rationalisation in terminology. 
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Likewise, there is no doubt that metamorphic textural terms are often 
genetically appropriate, but there must be occasions when that particular 
origin is doubtful. Anhydrite beds at depth sometimes grade into secondary 
gypsum beds at outcrop presumably, but not always certainly, as a result of 
weathering (Pettijohn, 1949, p. 359). Surely, unless one can be reasonably 
certain that no genetic implication is involved, it is unwise to use a meta- 
morphic-type description for those secondary gypsum beds. It is precisely to 
avoid problems of this type that the proposed system may be most useful and 
can be positively justified. 
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Handbuch der Physik. vol. 47, Geophysics 1: The Earth’s Body; |. S. 
Fiucee, editor. P. 659. Berlin, 1956 (J. Bartels, Group Ed. Springer Verlag. 
DM. 118).—Geophysics I is the first of two volumes dealing with geophysics 
in the new edition of the Encyclopedia of Physics. The volume consists of 16 
chapters: one each on rotation of the earth, radioactivity, telluric currents, 
interior of the earth, and forces in the crust; two on magnetism; two on 
gravity and figure of the earth; and seven on related topics of seismology, 
structure of the crust, and elastic properties. The book serves as an excellent 
summary of the fields of geophysics concerned with the solid earth. 

The volume opens with a clear and concise treatment of the rotation of 
the earth and the derivation of a constant time unit by Sir H. Spencer-Jones. 
Irregular variations in the rate of rotation are explained by Munk and 
Revelle’s theory of thermal convection in the core and electromagnetic coupling 
of the core with the mantle. 

The second chapter, the first of several on seismology. is by J. Coulomb 
(in French) and treats the theory and mechanical principles of seismographs. 
In the next chapter K. E. Bullen develops the basic equations of seismic wave 
transmission, discusses the cause and energy content of earthquakes and in- 
terprets the physical properties of materials in the deep interior of the earth. 
Next follows an outstanding chapter on surface waves by M. Ewing and F. 
Press. This chapter summarizes much new material that has never before been 
brought together. Geologists will be interested to note that Lg waves are in- 
terpreted as channel waves in a low-velocity channel the top of which is the 
surface of the earth, instead of requiring a rather special low-velocity channel 
inside the crust. The chapter also provides necessary background for a chapter 
on structure of the earth’s crust by the same two authors. In the latter chapter 
seismic refraction studies and gravity measurements are also brought to bear 
on the problems of crustal structure in the continents and oceans. The conti- 
nental crust is about 35 km. thick and the oceanic crust only 5 km. The 
favored theory on oceanic trenches associated with island arcs is depression of 
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the crust without thickening or development of a “tectogene.” Other chapters 
on seismology and related subjects are one by J. Coulomb on microseisms (in 
French), a very short chapter on seismic prospecting by Ewing and Press, 
and a chapter on measurement of elastic constants of rock. The latter, by 
H. Baule, Bochum, and E. Miiller (in German) reviews static, dynamic and 
seismic methods of measurement. 

In a chapter on gravity and isostasy, G. D. Garland describes absolute 
and relative measurements, their reduction, and their mathematical interpre- 
tation. An especially interesting section on the meaning of anomalies includes 
a description of the Bouguer anomaly on the Canadian shield, where the 
anomalies (with a relief of 90 milligals) follow structural trends and seem to 
indicate major lateral changes within the crust rather than simply variations 
in thickness of the crust. Correlation of negative anomalies with glaciated 
regions is also challenged. The subject of gravity and geodesy is completed in 
a chapter by K. Jung on the figure of the earth (in German). This topic 
covers nearby a hundred pages and includes a thorough, lucidly illustrated 
discussion of the geoid, the figures of reference and movements of the earth’s 
axis of rotation. 

The difficult and elusive subject of forces in the earth’s crust is developed 
hy A. E. Scheidegger, avoiding “obviously fantastic postulates.” Some readers 
may take exception to basic premises in this paper, for example the premise 
that mountain building occurred in only about ten diastrophisms “separated 
by long periods of acquiescence,” or the assumption of a two-layer crust fol- 
lowing Jeffries’ 1939 model, or especially, the claim that gravity data elimi- 
nate the possibility of normal faults along the African Rift Valleys. The latter 
premise is widely recognized as probably in error (see for instance, Vening 
Meinesz, F. A., 1950, Les “graben” Africains, résultat de compression ou de 
tension dans la croiite terrestre?: Institut Royal Colonial Belge, Bulletin des 
Séances, v. 21, p. 539-552). A comparatively long article on radioactivity and 
age of minerals by J. T. Wilson, R. D. Russell and R. M. Farquhar covers the 
subject historically as well as analytically, A chapter on the earth’s interior by 
J. A. Jacobs is concerned with seismic data, mass and moment of inertia, 
thermal and electrical effects, and the geomagnetic field as explained by the 
dynamo theory. Electrical currents in the earth are treated in an article by L. 
Cagniard (in French). The chapter includes general theory, electrical pros- 
pecting and relations between magnetic and electrical variations. 

Magnetism is treated in two chapters, both by S. K. Runcorn. The first 
is a concise summary of rock magnetism: mineralogy, physical theory, meas- 
urement, and variations from the Precambrian to the present. The other 
chapter summarizes the magnetism of the whole earth and theories of the 
earth’s magnetism. 

A useful German-English and English-German index completes the book. 

The volume is attractive, except for its high price, to specialists in geo- 
physics and to everyone interested in geophysical problems. 

GEORGE A. THOMPSON 
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Elasticity, Fracture and Flow, with engineering and geologic applications ; 
by J. C. Jaecer. P. viii, 152. London, 1956 (Methuen & Co., Ltd., $2.50). 

This splendid little book by the author of “An Introduction to Applied 
Mathematics” and co-author of “Conduction of Heat in Solids” is a concise 
mathematical treatment of the behavior of stressed solids. Professor Jaeger, 
in his own words, has “attempted to set out in as elementary form as possible 
the basic mathematics of elasticity, plasticity, viscosity and rheology, together 
with a discussion of the properties of the materials involved and the way in 
which they are idealized to form the basis for the mathematical theory”. 

The book treats three topics in approximately equal length: stress and 
strain, behavior of actual materials, and equations of motion and equilibrium. 
Stress and infinitesimal strain are presented first for two and then for three 
dimensions. Invariants, quadrics, Mohr’s representation of stress and strain 
then follow, with a careful elementary treatment of finite homogeneous strain. 
This section is extremely good although the notation is not best suited for 
extension to tensor analysis. 

The behavior of actual materials is described in an extraordinarily com- 
pact fashion which nevertheless achieves the stated purpose of the book. The 
equations of isotropic elasticity are followed by review of crystal elasticity, as 
an example of anisotropic elastic behavior. Treatment of finite strain under 
hydrostatic pressure next suggests to the reader what can be done with finite 
elasticity. Fracture of solids is treated according to Coulomb-Navier theory, 
with Anderson’s extension to faulting. The Mohr, Griffith and Becker theories 
of failure are summarized, although unfortunately the reader is not given 
references to the original work. The yielding of perfectly plastic solids accord- 
ing to the Tresca and Huber-Mises theories conclude this section, with a brief 
summary of rheological models which combine elastic, plastic and viscous 
behavior; 19 references are provided. 

The final third of the book develops the equations of motion and equili- 
brium of isotropic bodies, and uses the stress function to analyze a few prob- 
lems of plain stress and strain. Wave propagation, a few elementary relation- 
ships in seismic prospecting, the equations of motion of a viscous fluid, and 
the slip line theory of two-dimensional plastic flow of ideal materials conclude 
the book. 

Geologic applications briefly mentioned are Anderson’s theory of fault- 
ing, Becker’s finite strain theory, elastic properties of rocks, stresses at depth 
due to gravity, seismic prospecting and creep of rocks. Engineering applica- 
tions include, in addition, the tensile test, plastic torsion of cylinders, and the 
Rankine theory of earth pressures. 

This book should be useful to geologists who want a rather complete 
elementary survey of those mathematical tools presently available for analysis 
of rock deformation, or who want to review the assumptions one must make 
about rocks in order to treat them mathematically. For mathematical prepara- 
tion the reader should have facility with calculus through partial differentia- 
tion and vector calculus, and some acquaintance with three-dimensional 
analytical geometry and the complex variable. The book is nicely subdivided 
for classroom use; an elementary course in structural geology could use sec- 
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tions 1, 2, 3, 6, 7, 10, 21 and 22, which treat two-dimensional stress and 
strain, Mohr circle and strength theory, and general properties of materials. 
Remaining sections might be used in an advanced course in mechanics of rock 
deformation or crystal plasticity or as an introduction to seismology. Com- 
pactness of the book and the need for supplementary reading and problems 
mean that this book will find its greatest application in the hands of a teacher. 
Probably the greatest weakness of Professor Jaeger’s book is that it is in 
places forced to touch on the physics of deformation, but for lack of space 
cannot really survey this field. For example, the most recent work reported 
on rock or mineral deformation is from 1939, the rich field of crystal plastic- 
ity is totally unrepresented, and thermodynamic variables are not considered. 
Furthermore it is implied that with present knowledge we can idealize the 
flow of rock by a simple rheological model (p. 105). In this respect the title 
is somewhat misleading for there is regretably little in geology that can now 
be treated by the methods of this book without further experimental guidance. 
However, the book succeeds admirably in its expressed intent, and as a well- 
written synthesis of material now scattered through six or seven books on 
plasticity, elasticity, metallurgy and beam theory, it is highly recommended. 
W. F. BRACE 


Principles of Stratigraphy; by Cant O. Dunpar and Joun Ropcers. P. 
356; 123 figs., 21 tbls. New York, 1957 (John Wiley & Sons, Inc., $10.00) .— 
The present status of stratigraphic principles is ably and clearly described by 
Dunbar and Rodgers in Principles of Stratigraphy, a readable, scholarly sur- 
vey well-suited for use as a college textbook and equally valuable to the pro- 
fessional geologist as a reference work. It is comprehensive, up-to-date, and 
attractively presented, the outstanding book in the field. 

Eighteen chapters are divided into four sections as follows: Environments 
of deposition (sedimentary processes, non-marine environments, marine en- 
vironments, mixed environments), basic stratigraphic relations, (stratification, 
breaks in the record, facies and facies change), interpretation of lithotopes 
(sedimentary rock nomenclature, rudites, terrigenous arenites, lutites, redbeds, 
carbonate rocks, siliceous non-fragmental rocks), synthesis (the local section, 
correlation, the stratigraphic system, broad patterns in the distribution of 
sedimentary rocks). A combination of bibliography and author’s index, listing 
685 titles, and a subject index are included. 

The method chosen by the authors for dealing with inconsistent and 
multiple use of terms is logical. It should go far to resolve semantic conflicts 
that abound in stratigraphy. They believe that “what can be dangerous to 
scientific understanding and progress is not so much that a term is being used 
in two senses, annoying as that may be, as that the men so using it are un- 
aware of their difference.” They recognize that uniform usage, although 
desirable, cannot be forced. So, they describe differences and inconsistencies 
that exist and make their own recommendations. 

A few minor points invite comment. Microcrystalline texture in carbonate 
rocks is, in some cases, clearly related to fine dimensions of the original 
grains. However, there are well-authenticated cases (e.g., in the Guadalupe 
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Mountains area, Texas) of development of secondary microcrystalline texture 
in carbonate rocks by recrystallization of coarse clastic sediment. The clastic 
texture and fabric of the original rock are progressively obliterated by re- 
crystallization. 

The authors devote much space to recent studies on turbidity flows, but 
they do not make a clear distinction between them and submarine slides 
(termed submarine “landslides” ). The slides may cause turbidity flows, as has 
been postulated for certain Permian rocks in the Delaware Basin, but the 
deposits of the two processes are very different. Turbidity current deposits are 
characterized by graded bedding and, according to Kuenen, the coarsest frag- 
ments do not exceed 10 cm in diameter. It is doubtful, and certainly not 
demonstrated, that turbidity currents “are capable of transporting very large 
blocks.” 

The authors use the term biotope Tor the environment in which a com- 
munity of organisms lives or lived in the past. They do not distinguish, how- 
ever, between organic environments (habitats) and depositional environments. 
It seems to me that these are usually different and should not be confused. In 
ecological usage, the biotope is a map unit. It is synonymous and coextensive 
with niche or habitat, But many habitats are not closely related to conditions 
of sedimentation. By simple analogy one would suppose that a lithotope would 
be the depositional (and diagenetic) conditions under which a rock was 
formed, and it is so used by some writers. However, Dunbar and Rodgers 
define this term as the local rock record of the biotope. 

Biocoenosis, another ecological term, is used by the authors in the con- 
ventional sense for a local community or assemblage of living organisms. If 


this definition is to be maintained, however, it cannot very well be applied 


to fossil assemblages which bear little resemblance to living communities. 
Very few species of the living community, perhaps rarely more than a few 
percent, are capable of leaving a fossil record. Fossil assemblages are records 
by which biocoenoses may be inferred. 

A few points on coral reefs raised by the authors are debatable. They 
subscribe to a currently popular view that the broad sea-level flat, so char- 
acteristic of Pacific coral reefs, is cut mainly by solution. There is, however, 
an opposed view supported by much evidence that the reef flat is of complex 
origin in which sea-water solution is subordinate. According to this thesis 
planation near the inner margin of the reef is by organic erosion, aided by 
wave attack at the intertidal nip. Elsewhere, over most of the area of the reef, 
the intertidal platform is built upward to sea level by organisms. The surface 
represents the ceiling of organic accretion. This, essentially, was the conclu- 
sion of Darwin and Dana and is the view shared by a number of modern 
students. 

The Marquesas Islands near Key West cited by the authors as an atoll 
are not of coral reef origin. They form a pseudoatoll. 

These are relatively unimportant details chosen at random. They do not 
detract significantly from the general effect of objectivity and comprehensive 
treatment of Principles of Stratigraphy. The authors and publisher are to be 
congratulated for the outstanding and timely result. NORMAN D. NEWELL 
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